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The general problems were investigated by the grantee in an 
lytical framework. These are the location of transfer stations a 
routing of collection vehicles. Mathematical optimization models 
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algorithms. These models, particularly the ones for transfer fac 
site selection, may be used in choosing the economically optimal 
from among a large number of alternatives. This information serv 
nonquant if iable information in making the final choice. 

Several of the models were tested by the grantee using data 
from the City of Baltimore. Feasibility of transfer stations, th 
sirability of rail haul, and the cost of increased collection fre 
were investigated in the testing process. These applications ser 
demonstrate the final usefulness of the models. 
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ABSTRACT 


The application of operations research to the i 
solid waste collection systems is studied. Models an< 
for facility location and routing are discussed and e: 
Three different types of problem areas are viewed. Cl 
covers facility location, in particular the location < 
facilities within a large-scale system. Chapter III t 
finding optimal flow through given systems with added 
In particular, a multi-commodity truck assignment pro! 
a common vehicle fleet is used to carry several commo< 
between supply and demand points is extended and solvt 
Chapter IV is concerned with vehicle scheduling prohle 
routes are to be found for individual collection vehi< 
various tasks. In Chapter V, the analysis of an actui 
scale solid waste collection system, that of Baltimor< 
is carried out using some of the methods developed. I 
the analysis indicate that the use of such models in 1 
of large-scale public systems can supply a great deal 
information about, and insight into, the management ar 
of the system. 
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CHAPTER I. OPERATIONS RESEARCH AND SOLID WASTE COLLECTION 

reduction 

This thesis is concerned with the development of analytical 
hniques to aid in the management and planning of solid waste 
lection systems. Specifically, the techniques to be presented 
operations research models of routing and facility location 
blems that are an inherent part of any collection or delivery 
tem. The models are most applicable in the private sector of the 
nomy, since they seek to minimize a well-defined and quantifiable 
t subject to explicit constraints. Solid waste collection 
terns, however, are for the most part functions of the public 
tor rather than the private sector. There is a divergence be- 
en the objectives of a public and private system, with the public 
ective function being more vague and difficult to express formal- 
constraints on the public system, especially those of a 
itical or social nature, may be difficult to measure, and criteria 
effectiveness may not exist in units which are commensurate. 

All this does not mean that operations research models and 
Is should not be applied to public sector problems such as solid 
te collection system management. As with private sector models, 
lysis is performed as an aid to the decision maker and not as a 
laeement for him. Careful use, based on an appreciation of the 
tors that the models cannot consider explicitly, will allow 
ningful analysis to be carried out. The amount of analysis 


attempted in such large scale systems until recently has been 
limited by the sheer size and complexity of the problems. Now w: 
the availability of operations research techniques and the aid o; 
high speed computers, the potential to increase greatly the undei 
standing of such systems is apparent. It is hoped that this wil. 
be an essay both in the advantages to be gained by systematic 
investigation and the pitfalls involved in blindly following modi 
as a substitute for judgment and intuition. 

The purpose of this work is twofold. First, the routing an« 
facility location problems of a collection system will be analyz< 
and specific techniques suggested for approaching their solution 
Second, the application of these techniques to the investigation 
a large scale solid waste collection system as an indication of i 
questions amenable to analysis by these models will be presented, 

A. The Scope of the Solid Waste Problem 

The solid waste system is defined to include the generation 
waste at a source, its collection and transport and the disposal 
methods available. It is a large, expensive and vital on-going 
system with peculiarities and problems that are in urgent need o: 
analysis. The yearly cost of operating the system in the United 
States is reported by Black (1964) to be 3 billion dollars per yi 
This amount makes it one of the foremost public works expendituri 
ranking after education, highways, welfare, and fire and police 
among local government expenditures. The amount of waste genera 



the United States from all sources approaches 125 million tons 
: year. With increasing indications that these costs and volumes 
LI continue to rise drastically (Ludwig and Black, 1968), it seems 
it detailed analysis of the system would be in great demand. 

The main difficulty in analysis is determining some measure of 
fectiveness for the service provided. The evolution through the 
irs of solid waste management from a private to a public problem 
indicative that the service provided is a public good. The 
;ire community suffers from the viewpoint of public health and 
sthetics when even one of its members refuses to dispose properly 
his solid wastes. The cost of enforcing sanitation ordinances 
:hin a private system may far outstrip the cost of a publicly 
led system. A measure of efficiency of a public system might be 
> degree of service offered to the customer in terms of frequency 
collection, types of wastes removed, locations from which waste 
collected, and the general level of satisfaction shown by the 
isumers in order to provide the incentive for all to dispose 
>perly of waste material. 

Management to achieve some desired level of service at minimum 
it is the goal, then, of the system operator, and he must study 
5 various control alternatives available to improve the efficiency 
the system. Cost analyses by Ludwig and Black (1968) reveal that 
percent of the solid waste system cost is due to collection and 
Ly 15 percent to disposal. This does not necessarily mean that 
2 only approach to the problem is through improvements in 


collection. In spite of the fact that most analysts consider tfo 
waste loads are given and begin designing a system from that po: 
it would seem that collection costs could be decreased by elimin* 
wastes before they require collection and transport. From a tecl 
logic point of view, materials which are difficult or expensive i 
handle or dispose of might be replaced by those which are not. ' 
impetus for such a move could be either legislative, through laws 
preventing the use of objectionable material, or economic, throuj 
special taxes which would make using the material economically 
infeasible. Economic incentives might also change current waste 
material into production inputs. Education of the public might i 
reduce waste load produced. Technologic advance in the form of 
better in-home disposal units is a hope for the future even thoug 
present devices are less than adequate. Present day incineratorf 
create an air pollution hazard and still have an ash disposal 
problem. Garbage grinders simply trade a solid waste collection 
system for an underground liquid waste system, and disposal probl 
continue. 

Some ideas have already been advanced for improving collect^ 
Prominent among these is the building of transfer facilities witf 
a city at which the special purpose collection vehicles could dis 
charge their loads to special purpose transport vehicles and thus 
return to their collection jobs sooner. The suggested transport 
vehicles have ranged from large tractor trailers to barges and r£ 
road cars. Recommendations have also been heard for even more 


ialized collection vehicles such as small carts hauled in tanden 


jeep or motor scooter to a rendezvous with a larger transport 
cle, Zandi (1968) has proposed an extensive pipe line system 
moving solid wastes under pressure as an alternative to surface 
ection. All of these control methods require some amount of 
stigation and testing. Yet, the most important thing to realize 
t the solid waste system is that it is too big, complex, and 
1 to allow actual experimentation without great expense or the 
ntial of great chaos. Coupled with this are all the other 
lems involved with studying large scale public systems. A rele- 
data base is probably non-existent. The political implications 
ho controls the system and who pays for which service may be the 
persuasive argument for throwing out a scheme that might 
rwise seem quite efficient. Since large investment has already 
made in attempting to manage the system, the designer is deniec 
luxury of starting from the beginning and is often saddled with 
ing around the blunders of the past. To this morass the 
yst hopes to bring some order. 

Since collection makes up 85 percent of the cost of the solid 
e management system, it is to this phase that the attention of 
work will be directed. Successful analysis of the collection 
ess could result not only in means for making it run more 
clently, but also provide a means for measuring how non- 
ection alternatives which reduce waste load must be priced in 
r to make them better choices than the current approach to the 
lera. 


B. Literature Review of Operations Research in Solid Waste Col 


A search of the pertinent literature reveals that very few 
authors have concerned themselves with the application of opera 
research to the study of solid waste collection systems. The w< 
that does appear falls into two categories: broad scale attemp 
outlining the nature and interactions of the entire waste manag< 
system, in which collection is but one part, and narrower, more 
specific investigations of the collection and disposal operatioi 
suited to the detailed study of present operating systems. The 
complete of the broad scale studies is one still under way at tl 
Univecsity of California at Berkeley entitled "Comprehensive Sti 
of Solid Waste Management". In the first annual report of the t 
(Golueke and McGauhey, 1967) one of the stated objectives of th< 
study is "to explore the potential of operations research to he: 
in the definition and solution of the solid waste disposal prob! 
Basically, two problems were investigated. The first was to fii 
means for evaluating solutions for the total refuse disposal prc 
that would not only be applicable to a wide spectrum of alternaJ 
but would be sensitive to environmental and governmental constrs 
as well. This involved the development of conceptual models whi 
gave insight into the operation of the system. As a second prot 
the study undertook to build mathematical models for the concept 
models. These included a waste generation model, a waste collec 
treatment and disposal model, a regional economic model, and moc 



for population, public health aspects, land use and process te 
logy. The most interesting of these models from the viewpoint 
operations research techniques was the waste collection, treat 
and disposal model. The subgroup working with this problem lo 
at the means of determining how solid wastes should move from 
sources through treatment and processing to disposal sinks. T 
identified the problem as one of network analysis and noted th 
solution should be gained by a graph theoretic approach. The 
come of the development of specific techniques for solving the 
as formulated is reported by Anderson (1968), who was one of t 
study group. Anderson was interested in modeling how the flow 
be optimized through a given system which contained existing 
facilities and potential new facilities. He identified the pr 
as a study 9 f a trans-shipment network with some peculiarities 
to the nature of the processes involved. In considering a tot 
waste management model including intermediate treatment facili 
it is necessary to consider the changes in volumes and in the 
of the wastes which is not possible using ordinary network alg 
He suggests two solution techniques: if the number of disposa 
is small, a branch and bound model may be used to generate opt 
flow through the system. If not, he also presents an out-of-k 
algorithm for solution which includes an added feature of alio 
flow through a node to be proportioned in a fixed ratio. In b 
models all costs are linear, and the system under investigatio 
considers the gross shipment between points an.d not the routin 


individual vehicles among small collection tasks. Also, 
dealt with the problem of finding flows through a known 
configuration. To consider various combinations of new 
each scheme must be evaluated separately. No methodolog 
presented for a systematic search over facility location 
ties. The assumption of linear costs omits some of the 
aspects of the problem. Although nonlinear costs may be 
linear approximations, this is only valid if the nonline 
are convex, Facilities costs are commonly concave or qu 
(fixed charge). Thus linear approximation either will n 
successful or must be clone in such a manner that one los 
confidence in the meaning of the results of the solution 
Other studies defining the solid waste system have 
out by consulting organizations for governmental units, 
any have gone beyond the characterization of the system 
at simulation of the process. Such studies include repo 
Management Technology, Inc. (Anon. 1966b) for the State o 
and by Aerojet General (Anon. 1965) for the State of Cal 
The first work in attempting to understand the corap 
microstructure of the detailed collection process at the 
block and vehicle level has been by simulation models, 
and Wersan (1965) presented a model of a collection syst 
particular attention to queuing problems at the disposal 
using data from Winnetka, Illinois. Truitt, Liebman and 
1969) constructed a simulation model in which the effect 


proposed transfer station may be investigated. Their results 
based on data from Baltimore, Maryland. Quon, Tanaka and Wer 
(1969) continue their earlier simulation work with a model fa 
studying changes in work rules and collection policy. 

Other authors have dealt with different aspects of the c 
tion problem that are more amenable to analytical attack. Co 
and Martin (1967) presented a heuristic method based on dynam 
programming for aggregating small collection areas into work 
schedule assignments for crews and vehicles. Shelly (1968) p 
ted a fixed charge model for looking at the large scale probl 
transportation of wastes in which the variables include the a 
native location of transfer stations. He used for a solution 
technique a heuristic fixed charge algorithm developed by Wal 
(1968). Skelly's work is discussed in greater detail in Chap 
Wolfe and Zinn (1967) have presented a simple model for an ov 
crude evaluation of a large system as an example of systems a 
cation in public works problems. A small trial problem was s 
using linear programming. 

C. The Analysis of Solid Waste Collection 

The analyst involved with a solid waste collection syste 
concerned with answering a set of questions about the system 
he can better understand it. This understanding will lead to 
improved decision making. The sort of questions he might ask 
the following: 



1 . 


What are the goals of the system? What frequei 
collection and types of service should be offei 


system? How will changing the service affect c 

2. What types of vehicles should be used, and how 

3. How many personnel are needed and what should t' 
duties and work rules be? 

4. What route should be assigned to each vehicle? 
the city be divided into administrative subgroup 

5. Are there parameters of the system to which syet 
and variables are particularly sensitive? 

6. If there is additional money available for resea 
what aspect of the system should further study bi 

7. Should there be intermediate transfer stations fc 
deployment of wastes to more specialised transpot 
Where should they be located and what type of equ 
should they contain? 


^ ^ UQCU 


8. What type of transport vehicle would 

from a transfer station to final disposal? 

9 ' ° f diSp ° Sal ^-native should be c 

where should it be located? 

10, ““ W0UW b6 <>■* the system of new 

in in-house waste reduction? m new disposal 

How will the stochastic nature of waste genera 

the 8nalySlS? H0W Wl11 *• solution change as 
be served continues to grow and spread? 


12, What are the effects of political, social and economic 
constraints? How much should be spent on aesthetic 
factors? Is regional grouping a feasible alternative? 

Ideally, to answer all these questions, the analyst would 1 
to build a model that would be able to follow the collection prc 
in its minutest detail and be able to manipulate all of the myri 
parameters that could possibly affect the solution. However, it 
magnitude and complexity are such that to consider the system in 
entirety and encompass every possible detail is quite difficult, 
not impossible. For this reason, simplifying assumptions will fc 
made in model development which do not allow the detail mentlone 
but will allow some approximation of the problem. Basically, tv; 
different approaches are taken. If one is willing to neglect th 
problems of routing of the individual vehicles through their col 
tion tasks, the large scale problem of how material should move 
through the system and how transfer facility location may be chc 
can be approached optimally. This will be called the "tnacro-sca 
or "flow of materials" subproblem. If the location of facilitie 
and the collection areas assigned are known for a set of collect 
vehicles, the small scale problems of routing vehicles among 
individual collection tasks may be approached. This second sub- 
problem will be called the "micro-scale" problem. Extensive 
discussion and development of both of these problems will be pre 
sented. Application of the proposed models to the investigatior 
a solid waste collection system will be emphasized. 


The general outline of this thesis is as follows: 
is concerned with the location of transfer facilities. 1 
problem is treated as a large scale "flow of materials" \ 
is found to belong to a class of problems known in the li 
as location problems, warehousing problems and plant loci 
problems. Techniques for solving these problems are die< 
a special method is developed that is particularly adapt* 
so Lid waste collection problem. 

Chapter III also deals with the large scale flow < 
problem. In this case, the network structure of the sysi 
and optimal flow through it is found for more extensive < 
than were allowed in Chapter II. In particular, multicot 
using a common vehicle fleet is studied. Extension of a 
xnodity truck assignment problem by Szwarc (1967) is pres< 
computational experience is discussed. 

Chapter IV deals with the micro-scale routing prot 
are found to be represented by several classes of problei 
tions research. The most general problems discussed are 
vehicle scheduling problems and traveling salesman probli 
original algorithm for the m-salesman traveling salesman 
presented. Problems of routing in street networks are sJ 
cases to be Chinese postman problems, and this problem ii 
Chapter V deals with the analysis of an actual sy* 
of the techniques developed in the thesis with fairly cri 
Chapter VI presents a summary, and some conclusioi 
suggestions for further research. 



CHAPTER II. THE LOCATION OF TRANSFER FACILITIES 


A, Problem Description and Literature Review 

Determining the location of intermediate facilities where 
transfer of solid waste may take place from specialized collec 
vehicles to vehicles more suited for long-haul transportation : 
problem that has received considerable attention within the pa< 
years. The fundamental questions involve the desirability of 
transfer stations, their number, location, and capacity, and tl 
specific functions which they should perform. These decisions 
be viewed as a trade-off between the building of facilities anc 
cost of transportation. Two points are immediately clear. Fii 
this is a general problem of facilities location in both the pi 
and public sector rather than a problem specific to solid wastt 
and studies in the general area of facilities planning may be 
applicable. Second, assuming that facility costs and transpori 
costs can be defined, it may be possible to derive some mathemi 
formulation of the problem. 

A search of the operations research and economics literati 
reveals that a considerable interest has been generated and mat 
matical formulation for certain aspects of the problem has beei 
developed. The methods developed have been for the large scalf 
"flow of goods" features of the problem in which the concern hi 
been with selecting paths over which materials flow from colle< 
areas to facilities for transfer and disposal. The routing of 



Lcles within collection areas is usually ignored, as it would 
a complexity to the problem that would make it intractable, 

>ver, since the goal of model formulation is to obtain some 
jcture that may be solved and manipulated in order to gain better 
jrstanding of the system as a whole, such work is still of great 
>rtance. Care should be taken, though, in the use of the solu- 
is generated. An optimal answer to a subproblem may have little 
ling in the context of the larger problem, even though manipula- 
i of the subproblem to show the sensitivity of its optimal 
ition to changes in parameters and control methods may provide 
reat deal of insight and information for the analyst. 

The problem area as discussed in the general literature is 
by several names: the location-allocation problem, the plant 
ition problem, the central facilities problem and the warehousing 
?lem. To a large extent, all of these names imply the same 
jral problems, with differences depending on the assumptions made 
irder to facilitate modeling and solution techniques for the 
2ls. As noted, the general formulation is based on trading off 
Llity cost versus transportation cost. The problem is to find 
: system configuration that minimizes the sum of the two costs. 

Analysts who have worked on this problem have been faced with 
basic problems in definition and assumptions which have led to 
different structures for analytical approach. The first assump- 
\ is whether demand is spread, uniformly or otherwise, continuous - 
jcross the entire area to be served, or is located at discrete 


Lnts. The second assumption is whether potential facilities sites 
5 infinite in number or are limited to a discrete number of 
isible alternatives. If the assumption is made that both demand 
1 facilities locations are infinite the analysis procedure is 
Lte ill-defined and little analytical success has been obtained 
such problem solutions. However, once some form of limiting 
sumption is made of the solution space, considerable progress in 
Lution technique has been found. It should be pointed out that 
:h discrete assumptions are not severe from an engineering point 
view. In fact, within a city there are not an infinite number 
sites for solid waste facilities, but a finite number based on 
lilable land, zoning and concentrations of population and the 
ration of existing structures. Demand, too, does tend to cluster 
l an assumption dividing a region into areas with each represented 
its demand at some centroid is not unreasonable in most cases, 
j two approaches may be structurally distinguished by regarding 
5 problem with infinite solution space for demand and facility 
:ation as location on a plane, and the discrete problem as 
nation on a network. 

The problem names most associated with location on a network 
i the warehousing problem and the plant location problem. In the 
rehousing problem, warehouses or central facilities are built to 
:ve as supply points for satisfying demand in the region and the 
3 W of goods takes place from the facility to the demand area. In 
5 plant location problem, the central facility is a production 
rility which requires a flow of goods to it from the surrounding 


regions. From an analysis point of view, the problems 
and the problem names may be used interchangeably. Tl 
description of a warehouse problem may be converted tc 
of the plant location problem by changing the directic 
goods . 

The statement of the warehouse problem is: Giver 
demand areas for a certain product, each with a demand 
number of alternative sites where facilities may be bu 
these demands, determine where facilities should be bu 
demand areas are to be served by each facility. The o 
that the sum of the transportation cost and the amorti 
cost is minimized. 

The general mathematical formulation is: 
n m m 

Minimize; ^ + z F 1 (y 1 ) 

n 

subject to: i Xij = y ± i = l, 2 in 


ra 


2 X i1 “ D i 

1=1 J 3 

J-i. 

2, , . . . , n 

x i j ^ ® 

i=l, 

2, ..... m 


2, ..... n 

yi» o 

i-1. 

2, . . . . , m 


X ij 


where: 


- amount shipped from warehouse 
demand area j (a decision vari 



totaL amount shipped from warehouse i 
(a decision variable) 


= cost of shipping from i to j (dolla 


Ff(y^) » cost of building and operating warehous 
when y± is to be shipped from it (dolla 


D 


n 


m 


yi is to be shipped 
== demand at area j 
= number of demand areas 
« number of proposed warehouse sites 


If the functions F^(y i ) and are linear, the problei 

nes a simple transportation problem. Unfortunately, the funet 
l) is frequently nonlinear, and generally exhibits a large fix 
5tment for land, foundations, utilities, etc., before any amou 
>e manufactured or shipped. Once the facility is begun, the 
Inal cost of another unit of storage or manufacture may decrea 
icmies of scale). Such a function is termed concave and is no 
ible to linear programming. This is the feature which has dra' 
nost attention from researchers in this field, and the method ■ 
ting the concave cost function is the key element of most 
2 dures. 

Among the earliest work in the literature is a paper by Bau: 
folfe (1958), in which the authors approach the problem by usi 
*ear facility cost instead of a fixed charge. Then the problem 
?e formulated as a transportation problem. The resulting solu 
is then adjusted, using a heuristic procedure, to reflect the 
3 chafge aspects of the problem. This discussion will be 
:ed to five fairly recent approaches, all thought to have some 


application and validity. The papers by Feldman, Lehrer and Ray 
(1966) and by Kuehn and Hamburger (1963) represent two early forn 
lations of the problem which rely on heuristics to achieve good 
solutions, Balinski (1965) formulated the problem as an integer 
programming problem which is capable of solution for small cases. 
Efroymson and Ray (1966), Spielberg (1969) and the author of this 
thesis propose methods which achieve optimal solutions to larger 
mathematical problems as structured through branch and bound 
techniques. 

As noted before, Kuehn and Hamburger used a heuristic algori 
to find good solutions to the warehousing problem. They assumed 
that transportation costs are linear functions of the amount trar 
jorted and that facility costs are of the form, 

Fi( yi ) = if the facility exists 

=0 if It does not 

Thus F^(y^) consists of a fixed charge that is independent < 
:he storage and a linear cost which does depend on storage, if tt 
facility exists. 

The method starts with one facility and tries adding anothei 
facility to see if the total cost can be decreased. The general 
issumption is that the best N facilities contain the set of the t 
f-1 facilities. Termination of the procedure occurs when it appe 
:hat another facility cannot be added without increasing the tota 
:ost. 

One difficulty in working with heuristics is that there is v. 


r to tell how far the final answer obtained is from the true 
:imal solution. This difficulty tends to decrease the validity 
an analysis of sensitivity by obscuring the cause of a change in 
? solution which might occur when a parameter is altered. An 
jrovement could be due to the parameter change or to luck in the 
iputatlons. Thus, extreme care should be taken in using a 
iristic for sensitivity analysis. 

Feldman, Lehrer and Ray assumed a more general form for the 

:ility cost, that of a continuous concave function. They also 

turned that transport cost is linear in the amount shipped between 

.nts. Since F i (y i ) does not have a linear variable cost but a 

riable cost that is dependent on how many demand areas are being 

•ved by the facility, the problem of assigning demand areas to 

.sting facilities is not as easy as before. An approximation 

:hod is used to accomplish this assignment. Whereas Kuehn and 

iburger began with one facility, Feldman et al. started with all 

:ential facilities existing and dropped some out. Termination 

:urs when no further savings can be obtained by dropping facili- 

:s. Balinski (1964, 1965) developed a formulation for the plant 

nation problem that could be solved by integer programming tech- 

tues. In order to do this, F^(y^) is specified as a single fixed 

irge, F^, if the facility exists. Capacity constraints are not 

jcified for the facilities. The formulation is as follows; 

n m m 

l l c t . x tj + 1 y. 

j-1 i-1 j 1-1 


ilmize : 


( 2 - 6 ) 


subject to: 
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1 - i 

j — 1 , 2, . . . . , n 


0 $ x £ . i 1 

*-l. 

2 , . • . • , m 

j=l, 

2, . . . . , n 

yi = (0, l) 

i-l. 

2 m 


where: = t^ * 

t^j = the unit transportation cost from faeili 
to customer j 

Dj = demand at j 

x^j = the fraction of Dj supplied from facilit 
F^ a the fixed charge associated with facilit 
n = the number of customers 
m = the number of proposed facilities 

Equation (2-7) specifies that each customer must have a 
demands satisfied and inequality (2-8) requires that demand 
be satisfied by a facility unless it exists. Note that even 
small problem the above formulation has a formidable number 
stralnts and variables. For an N-customer, M-plant problem 
are NxM + N+ M constraints and N x M + M variables. Thus 
the case M a 10 and N=30, the problem to be solved has 340 con 
and 310 variables. 

It would seem, at first glance, that if the integer req 
(2-9) were replaced by the simple constraint £ 1, the sol 
the resulting linear programming problem would be integer, 
a counter example is reported and the means of solution firs 


impted is an integer programming cutting plane technique of 
>ry (1958). This is shown to have a very slow convergence even 
small problems, and is abandoned for larger problems. Balinski 
jests another scheme for solution which involves partitioning 
problem using a method of Benders (1962). This requires the 
:essive solution of more complex integer programming problems. 

•oof is supplied that the method will converge to the optimal 
ttion in a finite number of steps, but it appears that computation 
s will be quite large for this type of solution as well. 

Efroymson and Ray (1966) presented a method of solution to the 
• problem using a branch and bound scheme that has been tested 
mtationally with good results. Their formulation of the 
)lem is : 

m n m 

’mize: £ E c ii x ij + E F i ^i (2-10) 

i-1 j=l J i“l 

m 

ject to: E Xi, = 1 j=l, 2 n (2-11) 

i a l 3 


re: 


n 

l x i j ^ n 1 y i i-1, 2, m (2-12) 

K ij* F i ® (°* 1) (2-13) 

c. . =» cost of supplying the entire demand area j 
from warehouse i 

F^ = fixed charge for establishing warehouse i 
y^ » 1 if the i'th warehouse is built, 0 if not 


X. .= 1 if the demand area j is served b 
i, 0 if not 

n.. = maximum number of demand areas tha 

assigned to facility i 

m = number of possible warehouse sites 

n = number of demand areas 

An examination of the constraint set shows that Bali 
explicit constraints on the value for each allocation to 
(2-8) have been collapsed into a smaller set of constrain 
maximum amount that may be assigned to a facility (2-12). 
straint specifies that if the i'th facility does not exis 
no assignments may be made to it. However, if it does ex 
raents may be made to it up to the value of n^. The value 
becomes an important point in the development of the algc 
solution method involves an implicit search of all the pc 
feasible solutions. Suppose all the feasible solutions v 
sented as branches of a tree with each branch exhibiting 
value, zero or one, for each variable. In order to avoic 
the entire tree to find the optimal solution, means for € 
the lower bound on the solution by continuing down the t: 
available. If at any point on a branch the estimate of t 
bound is greater than a known feasible solution,, that brs 
dominated and may be abandoned. The speed with which sue 
procedure progresses to an optimal solution depends on th 
used to obtain lower bounds and the rules used for deterir 
which branch to investigate next. 


In establishing their branch and bound scheme to solve the 
problem, Efroymson and Ray had to specify a means for finding a 
lower bound at each node on a branch. They specified that this 
could be accomplished by solving the problem denoted by the obj 
tive function (2-10) plus constraints (2-11) and (2-12), while 
ignoring the integer requirement (2-13). If is a capacity c 
straint such that it is limiting, i.e. areas that might assign 
are excluded because there is not enough capacity at i, the pro 
to be solved at each node requires the solution of a transporta 
problem. However, if nj_ is greater than or equal to all the ar 
that could possibly assign to 1, the lower bound problem at eac 
node may be solved by a very simple iterative procedure that i 
very fast. Simply assign the area to the open facility such th 
transportation cost plus expansion cost at the site is minimize 
In order to speed their algorithm, the authors abandoned the 
requirement for a capacity constraint on facilities by setting 
at a number greater than or equal to the total number of demand 
areas that could possibly assign to that facility. Solution of 
problems with 50 possible warehouse locations and 200 demand ar 
in 10 minutes on the IBM 7094 is reported. 

Spielberg (1969), in a more recent paper, discussed largel 
the same problem. His algorithm for solution is also branch an 
bound, but contains many added features which shorten computati 
time. A great deal of computational experience is reported wit 
problems as large as 100 facilities and 150 demand areas. The 


author also suggested extensions of the problem, but reported no 
computational experience for these. 

Looking at the special problem of the location of transfer 
facilities for solid wastes, it is noted that there is an added 
dimension beyond that of the warehousing formulation. The facilj 
to be located are intermediate points within a trans-shipment net 
work between the sources and sinks of flow. Thus the transportat 
cost involves two elements, transportation from source to inter- 
mediate point, and from intermediate point to the sink. There ai 
possible facility costs at the intermediate point as well as 
variable costs at the sources and sinks. 

One of the weaknesses of all the previous work in warehousii 
problems is that there may be no capacity restrictions on the am< 
of flow through a central facility when in fact the cost estimati 
and the physical reality of the problem require that there be su< 
a constraint. The problem of locating an intermediate facility \ 
capacity constraints on flow will henceforth be called the capac: 
tated trans-shipment facility location problem. The only previoi 
work on this problem is that of Skelly (1968). In addition to 
recognizing the above features of the problem, Skelly was also 
interested in time-staging of construction of facilities and disj 
points. He formulated the problem with a series of linear const! 
which define the trans-shipment and time requirements. This lea\ 
a large fixed charge integer programming problem. 

Skelly chose as his solution method an algorithm by Walker ( 



r the fixed charge problem. This algorithm is simplex-based, i.e 
requires that a basis the size of the constraint set be estab- 
shed and manipulated. Furthermore, it is a heuristic and may not 
eld a global optimum. The fact that the Walker algorithm is 
mplex-based leads to some problems in computational efficiency, 
ans-shipment problems belong to a class of network flow problems 
described by Ford and Fulkerson (1962), which may be written in 
near programming form but because of the sparseness and special 
ructure of the constraint matrix, can be solved much more quickly 
d efficiently by primal-dual labeling procedures. Because of the 
mplex approach, the size of the problem that Skelly may consider 
severely limited. Furthermore, since the value of such models 
in the repeated solution using different parameters to show 
nsitivity, the use of a heuristic algorithm may be self-defeating 


B, The Capacitated Trans-Shipment Facility Location Problem 


The remainder of this chapter is devoted to a new solution 
method for the capacitated trans-shipment facility location prol 
where there are fixed charges associated with the intermediate 
facilities. The algorithm guarantees an optimal solution and w: 
solve fairly large problems within moderate computation time. ' 
problem to be solved is stated as follows: 

There is a set, K, of sources for a product with an amount 
at each source. In addition, there is a set of sinks, J, for t\ 
product, each with an upper and lower bound on demand of D^ u am 
A set of proposed intermediate facility sites, I, has been suggi 
as trans-shipment points between the sources and sinks. Each p: 
posed facility has a fixed charge, F^, a variable unit cost, 
which is a linear function of the amount shipped through the 
facility, and a capacity, Q^. The problem is to find which fac: 
should be built and which sources and sinks each facility servei 
that the total cost of facilities and trans-shipment is minimiz* 

In the context of a solid waste collection system, the floi 
through the system is solid waste material with the sources rep: 
ting the locations of waste generation and the sinks represent!: 
the disposal methods and locations. A lower bound on the araoun 
solid waste reaching a sink is motivated by the fact that for si 
types of present disposal processes such as incineration, a min 
throughput is necessary to keep the process operating. Future 



development In waste disposal technology may also lead to other 
processes which require continuous operation. The intermediate 
facilities are transfer stations where transfer may take place fr 
smaller collection vehicles to some other means of transport to t 
disposal points (sinks). The means of transport from the transfe 
station may be by vehicles such as large trucks or trains, or by 
non-vehicle methods such as pipelines. In addition to simple 
transfer, the activity which takes place at a transfer station ma 
include some form of processing such as sorting, compacting and/o 
incineration. 

In mathematical form this problem becomes the following! 
Problem One 

Minimize : £ F y £ + £ Z + Z Z c kl x kl (2- 

i=l 1=1 j=l i=l k=l 

subject to the constraints: 

^ x ki ^ S k kal) 2 * P 


n * P -ifif 

l x, , ~ I x ki i=l, 2, m (2' 

J-l J k-1 

P 

£ x ki ^ iml > 1 ' ra ( 2 

k«l 1 

u m * 1 

D >j * x A1 >/ D, J-l, 2, n (2 

J i=l J J 


where : 


* Vr* 

x ij * x ki are non-negative integers 


y± = (o, i) 


y± = 1 if the i'th facility is built 
= 0 otherwise 

x ij “ flow of material from facility l to sink 

irk 

X ki = flow of material from source k to 
intermediate point i 

* 

C ij “ c ij + r j * unit cost associated with a ti 
of material from facility i tc 
(dollars per unit) 

C lj = ““it shipping cost from facility i to sir 
(dollars per unit) 

r j = variable cost associated with using 
(dollars per unit) 

* 

C ki = c ki + c k + V i B uni t cost associated with 

transfer of material froi 
source k to facility £ 
(dollars per unit) 

C ki = “ nlt shl PPing cost from source k to facil* 
(dollars per unit) 

C k “ unit variable cost associated with using £ 
(dollars per unit) 

V i * unit variable cost associated with using 
facility i (dollars per unit) 

F i » fixed charge for establishing facility i 
(dollars) 

s k * amount supplied at source k 

Dj = upper bound on amount demanded at sink j 

D j X “ lower bound on amount demanded at sink j 

Qji * capacity of the i'th facility 



m 


number of proposed facility sites 


n = number of demand areas 

p = number of supply points 

Inequality (2-15) expresses the requirement that flow from the 
:e cannot exceed the supply of material there. Equation (2-16) 
2 S the conservation of flow requirement that the flow entering 
L*th facility must be equal to the flow leaving it. In 
lality (2-17) the fixed charge nature of facility location is 
sssed. If the i'th facility does not exist, = 0 and no flow 
rake place through it. If the i’th facility does exist, = 1 
flow up to may pass through it. Inequality (2-18) specifies 
flow to sink j must be between the upper and lower bound on 
rapacity of the sink. 

There are several additional constraints or side conditions 
might be of interest in the above formulation. A budget con- 
Lnt on the number or cost of the facilities to be built, or a 
il exclusivity constraint which allows a particular facility to 
I only if another does not, might be included. As noted with 
5 (disposal points) a lower bound requirement different from 
on the flow through a facility or a particular arc due to 
lical requirements might also be imposed. The side conditions 
be discussed in detail later in this chapter. 


G. Method of Solution 


It is desired to find a solution to Problem One that will gl 


values to the variables x^ ^ , x^ and so that the constraint se 


is satisfied and the objective function is minimized. One means 


doing this would be to set up the problem as an integer linear pi 


gramming problem and attempt to solve it by a cutting plane methc 


such as that of Gomory (1958). However, for an M-facility, N- 


deraand-area, P-supply-point problem there are 


2(M + N + P) constraints 


MxN+Mx P + M variables which are 


required to be integer. 

Thus for a problem of M=10, N=10, P=10, the integer program 
would have 60 constraints and 210 variables. Solution of the 


corresponding linear programming problem with the integer 
constraints relaxed will not necessarily yield integer solutions, 
However, there are some linear programming problems which h« 
special structure that may be related to a network, which may be 
solved optimally by procedures other than the simplex method of 
linear programming in a much shorter time and yield integer solu* 
tions. An example is the trans-shipment problem which represents 
supply-demand relationship where demands for a product at some 
points are satisfied by shipments from sources of the product 
through intermediate points. The resulting linear programming pi 
lem may be shown to have the form of a network, and finding the 


iraum cost maximum flow through the network gives the solution 
the trans-shipment problem. 

Consider Problem One. It has some elements of network struc- 
e, but it may not be represented as a network problem, 
qualities (2-15), (2-16) and (2-18) exhibit an incidence matrix 
ch by themselves would form a trans-shipment problem and would 
Id integer solutions. The addition of inequality (2-17) which 
a coefficient of one variable that may be different from zero 
one, destroys this property. However, a related but similar 
blem may be written in a form that will allow it to be solved 
a network algorithm. Consider the problem of determining the 
w through the system from supply source through intermediate 
ilities to demand areas if the facilities to be built are known, 
s implies that the values for the y^'s (the variable which 
ermines which facilities are built) are known and given. Let S 
al the set of indices representing the values of y^l. Then 
problem may be written in the following form: 
blem One (a) 

m n m p 

imize: I £ c-.x^ + z i c. + E F. (2-19) 
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where the variables are defined in the same manner as for P 
One, and 

S « the set of indices for which y^l. 


This problem has the form of a network problem as desc 

Ford and Fulkerson (1962), and can be solved by finding the 

cost maximum flow in the network. The method for finding t 

cost maximum flow through the network is the out -of -kilter 

1 / 

which was developed by Fulkerson (1961). One of the comput 
for the IBM 7000 and 360 series computers is in Fortran IV 

2 

and is available through the IBM Share Distribution System. 

It should be noted that without the upper and lower bo 
demand at the sink expressed in inequality (2-23), the abov 


- Those interested in following the development and proofs 
out-of-kilter algorithm are directed to Ford and Fulkerso 
or Fulkerson (1961). A good description of a more primar 
is found in Durbin and Kroenke (1967). A discussion of n 
techniques is found in Clasen (1968). 

— ^IBM Share Distribution #3536, "Out-of -Kilter Network or T 
tation Problem Solver", IBM, White Plains, New York. 


would be a simpler trans-shipment problem rather than an out-of 
kilter problem. However, the existence of the need for such 
constraints for this and other aspects of the problem means tha' 
an out-of-kilter problem rather than a simpler trans-shipment 
problem will be used for the basis for algorithmic development. 
The advantages of being able to formulate this subproblem as a 1 
work problem solvable by a network algorithm such as the out-of 
kilter are twofold. First, netwoi'k algorithms are very fast whi 
compared to simplex based linear programming techniques. Secon< 
the out-of-kilter algorithm may be started with any existing 
solution after constraints or costs have been changed. This re 
start capability will be of considerable interest in the soluti< 
scheme to be developed. Because of the network structure of th' 
problem, both the network algorithm and a simplex algorithm wil 
always yield integer solutions. 

The network structure of a location problem has been recogi 
by other authors. Efroymson and Ray (1966), and Spielberg (196 
while working on simpler location problems noted this structure 
took advantage of it in a branch and bound scheme. In their ca: 
the resulting structure was that of a transportation problem an 
the means of finding a lower bound on a branch might involve th 
solution of such a problem at each node of the branch. However 
the simplifying assumption is made that facilities do not have 
upper bound on capacity, the transportation problem decomposes 
a simple iterative scheme which is quite fast. Both authors 


attempted to save time by ignoring capacity constraints and usinj 
the simple iterative scheme. This limited their work to transpoi 
tation problems. 

The technique to be proposed below arose from the need to 
include both a trans-shipment network structure, and capacity coi 
straints on facilities in the formulation of transfer station 
location problems. It was clear from the outset that this more 
general formulation would require greater computation time, but 
it was anticipated that the increased cost of finding a solution 
would be more than offset by the greater generality of the soluti 
Much of the success of the algorithm is based on the ability to 
restart the out-of-kilter algorithm at each stage, so that littl* 
computational effort was necessary to calculate successive lower 
bounds. Trans-shipment and transportation algorithms, with sligf 
manipulation of the solution, also exhibit this restart capabilit 



A Sketch of the Algorithm 


It la desired to solve Problem One, which is repeated here for 
ivenience. 

ilem One - The Capacitated Trans-shipment Facility Location Problei 
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This problem has two elements that make it difficult to solve, 
•st, it is a fixed charge problem with the fixed charge variable, 
appearing in inequality (2-27) with a coefficient which may 
different from plus or minus one. Second, an integer solution 




is required. However, as has been previously noted, if the 
of the fixed charge variables, y^, are known, then the resu 

•f 

problem of finding the values of the flow variables, the x^ 
the becomes a form of trans-shipment problem whose s 

may be found by using the out-of-kil ter algorithm to find t 
minimum cost maximum flow of the corresponding network. In 
to solve Problem One, some approximations will be made. 

The flow of materials through the i'th facility is giv 


P 

Z 

k=l 


** 

x ki 


and the cost function of the i'th facility is 



if the i'th facility exi 
i.e. if y^-l 


0 if the i'th facility doe 

exist, i.e. if y^=Q 


The slope of this cost function is shown in Figure 2-1, alo 
a linear approximation of the cost. The approximation is a 
in the following manner. Convert the fixed charge, F^, int 
charge by dividing by the capacity of the facility Q^. The 
approximate cost function is 



and it underestimates the true cost function given in (2-29 
at two points where they are equivalent. These are when th 



FLOW THROUGH FACILITY 


.gure 2-1. The Fixed Charge Cost Function of Problem One and 
Its Approximation in Problem Two. 



and the use of bounds for eliminating areas of search, th< 
methods do differ somewhat in the search procedure. Impl: 
enumeration implies that a branching scheme will be set u] 
all possible feasible solutions will be enumerated either 
or explicitly. Bounds are used to avoid searching branch) 
are obviously dominated by existing feasible solutions. 1 
such techniques are Glover (1965), Balas (1965), Balas (l 1 
Geoffrion (1966, 1967). Branch and bound implies that th< 
feasible solutions is divided into mutually exclusive set) 
lower bounds estimated until an optimal solution has been 
and verified. In implicit enumeration, the goal is to eni 
all solutions while in branch and bound the goal is to aul 
into subsets until lower bound estimates show that an opt 
solution has been found. Examples of branch and bound tei 
are the work of Land and Doig (1960) and Little (1963). j 
of implicit enumeration in the field of environmental eng 
is a solution to a water pollution control problem by Liel 
Marks (1968). 

The method of solution chosen to solve Problem One i 
and bound. The basic reasoning behind the technique is t 1 
possible to solve a highly structured problem by solving 
structured surrogates. A feasible solution to the surrog. 
is always a feasible solution to the highly structured pr 
approximated as well. However, the value of the objectiv 
of the surrogate is not necessarily the value of the obje 


iction of the problem it approximates with the same feasible 
Lution. Because of the linear approximation used, the cost of 
i surrogate is lower than or equal to the cost of solving the 
Lginal problem. Define an exact solution as an optimal, feasible 
Lution to the surrogate where the value of the objective function 
: both the surrogate problem and the original problem are the 
ne, A lower bound solution represents the case when the cost of 
3 surrogate problem is lower than the cost of the original 
>blem. The surrogate problems are generated in such a manner 
it the solution space for the original problem is subdivided into 
:ually exclusive collectively exhaustive sets. The solution of 
i problem representing each set is a lower bound for that set. 
i optimal solution is found when the smallest of the lower bounds 
: all the sets is an exact solution to the original problem. 

Consider solving Problem One. The less constrained Problem Two 
LI be solved as a surrogate for it. A solution is generated by 
3 out-of-ki Iter algorithm and the cost and exactness of the solu- 
>n to Problem One are investigated. If, for each facility the 
>w through it is exactly zero or the capacity of the facility, an 
ict solution to Problem One has been found as well. Since this 
Lution represents the only lower bound estimate of Problem One, 
i it is feasible to Problem One, the optimal solution of the 
Lginal problem has been found as well. The algorithm terminates, 
rever, suppose that for some facility, i , the flow is greater 
in zero but less than capacity. Choose this facility to branch on 


In all feasible solutions to Problem One only two situat 
occur regarding i*. Either and the *• facility GX 
yi *=Q and the i* facility does not exist. Construct two 
problems: Problem Two (a), (i*)> is Problem Two plus the 
yi *=l. In order to set up this problem, the linear appr 
associated with i ^ i is replaced by its true 
is a constant, F,. , plus a variable cost of . No cha 
in the constraint matrix since it has already been assutr 
y^=l in the original formulation of Problem Two. Since 
constant, it is simply added to the cost of the solutior 
Problem Two (a), (i*), found by solving the resulting net 
the out-of -kilter algorithm. Since Problem Two (a), (i ) 
constrained than Problem Two, its cost must be greater t 
equal to the cost of Problem Two. Note that the flow t\ 
facility may take on any value between zero and capacity 
correct cost will be assigned, since the fixed charge fc 
facility is included in the cost, and the variable charg 
function of flow through the facility. Now construct Pi 
i , which is Problem Two plus the constraint “0. Thi 
handled in two ways. Either set the variable cost at f£ 
infinity or set the capacity at facility i to zero. B< 
dures will produce the result of blocking flow through 1 
The solution to the two problems represents lower boundf 
mutually exclusive sets of feasible solutions to the or: 
The lower bound with the lowest objective function is ei 


exactness to Problem One. If it is exact, it is optimal. If 
a new variable, i , is cho9en for branching, and two new prob 

^ iX. 

are set up. One problem has =1, =1, and the other has 

and y^ =0. The cost of these solutions, along with the remai 

lower bound from the first branch, are searched for a minimum 

more, and the exactness check is made. The procedure is conti 

until the optimal solution is found. The procedure must term! 

since the process of adding constraints continues only until t 

enumeration of all possible feasible solutions has been found. 

To present the subproblem solution in a general form, del 

S = the set of all indices representing facilities 

- the set of facilities constrained to be built 
(i.e. y^l) 

S° = the set of facilities constrained not to be bud 
(i.e. y^O) 

S ~ the set of facilities not yet constrained 
S* = S X \J S° u s 

s x n s°n s - 4> ; 

the sets are mutually exclusive and totally exhaust i-v 
Then a general subproblem may be written as: 

Minimize : 
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and 

that is , 


subject to the constraints: 
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Xjj, x kl are non-negative integers 


and all symbols are as previously defined. 


The ease of keeping track of the various solution 
puter is aided by the use of associative lists which a 
to be kept in order of increasing cost, and makes loca 
lowest cost solution relatively simple. 

The choice of a variable on which to branch can m 
deal of difference in the speed of finding the solutio 
used in the procedure for solving Problem One involves 
the facilities which have flow greater than zero but 1 
capacity of the facility. Denote the fraction of flow 
facility i as (defined as the flow divided by the c 
Select for branching that facility whose value of i 
to 0.5, which is the rule that Land and Doig (1960) su 


A Detailed Algorithm for Solving the Capacitated Trans-Shipmen 
Facilities Location Model 


The following data are given: 

m = the number of proposed facility sites 

n = the number of demand areas 

p “ the number of sources 

For each facility site, i, the following are known: 

= the fixed cost of establishing the i'th facility 

V ± - the unit variable cost of operating the 
established facility 

- the capacity of the facility 

For each demand area, j, the following are known: 

Dj U ® the upper bound on quantity demanded at 
the j * th demand area 

D. ^ the lower bound on quantity demanded at 
the j 1 th demand area 

rj - unit variable cost of satisfying demand at j 
For each supply point, k, the following are known: 

S^ u = the upper bound on quantity supplied at point k 
= the lower bound on quantity supplied at point k 
tj = unit variable cost of providing supply from point k 
In addition, the m x n matrix c^j , and the p x m matrix ar 
assumed known where 

c.. « the unit cost of transportation from facility i 
to demand area j 

» the unit cost of transportation from source k 
to facility i 



2. Initialization: 


Define S - the set of all indices of facilities 
have been fixed at zero or one 

S = the set of all indices not yet set 

S* » the set of all indices = (1, 2, 3, , , 

Then S U S = S* S f) S = (the empt 

Initialize S = 

S « S* » (1, 2, 3, m) « all indi 

seating 

3. Construct a graph in the following manner: There wi 
sets of nodes (I, J, K) which make up the graph, plu 

y 

node S#. The graph will be in circulation form so 
source = sink = S#. 

Define I » (set of nodes which represent the faci 
alternatives) 

For each of the m sites, a capacitated node is 
This is accomplished by representing the site 
with a directed arc joining them (see Figure 2 
upper bound capacity of the arc is the capacit 
the facility it represents. The lower bound c 
zero. The unit cost on the arc is Z^, which w 
defined later. 

Define J = (set of nodes which represent the dema 
Define K = (set of nodes which represent the supp 

3/ 

“"Circulation form means an arc exists from the sink 
the source. 



ttica uumiewLiug uit: uuuus; 


From S# to set K : Define a directed arc from S # to each 
node in K. The unit cost on the arc is t^. The upper boun< 
on capacity of the arc is S^ 11 and the lower bound is S^, 
Assume the upper bound is equal to the lower bound. 

From set K to set I : Define a directed arc from each node 
K to the receiving end of each capacitated node in I. The 
cost of the arc is c^ and the upper and lower bound capaci 
of the arc may bo set as desired. Assume the upper bound i: 
infinity and the lower bound is zero. 

From set I to set J : Define a directed arc from the trans- 
mitting end of each capacitated node in I to each node in J 
The unit cost of the arc is c ^j> anc * the upper and lower 
bound capacity may be set as desired. Assume the upper bon 
is infinity and the lower bound is zero. 

From set J to S# ; Define a directed arc from each node in 
to S #. The unit cost on the arc is r ^ , and the upper bound 
and lower bound capacities are respectively D^ u and Dji, 

Figure 2-3 shows a graph representation for m=3 , n“3 , k=2. 

Establish the unit cost, for the capacitated nodes in set 
Define « a decision variable 

Define y. = 1 if the i'th variable is in the set S and if 
1 the i ! th facility is to be built 


Define y. = 0 if the i'th variable is in the set S 
1 the i'th facility is not to be built 
if the i'th variable is in the set S 

Define Z i = 0 if y^»l and its 

Define - c© if y^O and i*S 

F, 

Define Z. = — =• if Vi*0 and IfS 

1 Qi 1 

5. Find the minimum cost maximum flow through the graph 
out-of-kilter algorithm. Define this cost as CGRAPH 
is the first time through this step and the attempt 
the graph using the out-of-kilter method indicates n 
solution, no feasible solution to the problem exists 
If not, continue. 

6. Define the total cost of the solution, COST, as 

COST * CGRAPH + E F,y< 
iCS 

7. Examine the solution to see if it is exact: 

Define X^ as the flow through the capacitated nod 
representing the i'th facility 

The solution is exact if for all ift£ 

= 0 or X 1 - 

If the solution is not exact, find the index of c 
members of the set S for which 0 < < Q^, as th 

candidate for branching. Select that member for 


nearest to 0.5. 



re information about the solution. File on a list ranked on 
cost of the solution information about the exactness, the 
bers of the S set and the values to which they are set, and 
next candidate for branching. 

ove the first problem on the list, i.e. the lowest cost 
ution found so far. If the solution is exact, it is 
imal. TERMINATE. 

the solution is inexact, set up subproblems (a) and (b): 

(10a). For subproblem (a): 

Set y^=l for the index stored as the candidate for 

branching and switch that index from S to S 

Then do steps 4 through 8 and return to do sub- 
problem (b) in (10b) 

(10b). For subproblem (b): 

Set y^=0 for the same index as in (a) and switch 
that index from £ to S 

Then do steps 4 through 8. 


Go to step 9. 
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F. Computational Results Using the Algorithm 


The algorithm as described was coded in ASA Fortr 
executed on an IBM 7094 computer. The out -of-ki Iter alg 
to solve the lower bound problem at each branch was a Sh 
Distribution Program written by Clasen; the list process 
routine package used in keeping track of the various sol 
written by Bellmore (1966). In Chapter V results based 
of data from example problems in solid waste collection 
General impressions of this type of algorithm gained by 
with randomly generated data are as follows: The usual 
branch and bound algorithms is that they are data depend 
same size problem with different data can lead to very d 
running times. In the case of the algorithm tested, the 
that ran the fastest for a given size were those where t 
costs were sufficiently high so that only a few faciliti 
built. From computational experience with such algorith 
been- found that problems where the optimal solutions are 
either a few or all of the facilities, generally run the 
Problems where there is a delicate trade-off between fac 
transportation costs, so that several facilities are bul 
take longer to derive optimal solutions. Further, the r 
out-of-kilter problems solved in obtaining the optimal e 
definitely exhibited the value of the algorithm's restat 
As the number of out-of -kilter solutions increased, the 


solution decreased. 


A possible limiting factor on the use of this algorithm is 
the amount of storage space needed within the computer for a 
lem is quite large. 

Let F = the number of proposed facilities 

D = the number of demand areas 
S = the number of supply points 
the number of arcs and nodes in the resulting graph is: 
Number of arcs =* D(1 + F (1 + S)) 

Number of nodes « D + 2F + S 

naximum size problem that could be run on the 7094 with its 
38 word memory was 300 nodes and 1500 arcs. The number of 
Lble alternatives is dependent on the number of facilities 
l investigated. However, comparing two solutions with the 
number of facilities but different numbers of supply points 
show that there is a significant difference in time to solve 
3ut-of-kilter algorithm due to its larger number of arcs. The 


itational results are shown in Table 2-1. 



Table 2-1. 


Results of Trial Runs of the Algorithm with Randomly Generated ] 


F 

D 

S 

No. of Al- 
ternatives 

No. of 
Arcs 

No. of 
Nodes 

No. of 
OKA's* ** 

70! 

Tii 

(mint 

6 

12 

2 

64 

104 

27 

8 

0 . 

12 

10 

2 

4,096 

194 

35 

8 

o.: 

* 10 

30 

3 

1,024 

382 

73 

8 

o.: 

* 10 

30 

2 

1,024 

382 

53 

8 

0.! 

15 

20 

3 

32,000 

388 

54 

3 

o.: 

30 

20 

3 

10 10 

733 

84 

15 

1.1 

10 

20 

3 

1,024 

293 

44 

8 

0 . 

20 

20 

3 

10 6 

503 

64 

6 

0 , 

30 

25 

2 

10 10 

866 

88 

30 

1 . 

+ 35 

10 

3 

10 12 

478 

84 

30 

1.1 

+ 35 

10 

3 

io u 

478 

84 

11 

O.i 

+ 35 

10 

3 

10 12 

478 

84 

6 

0 .. 

50 

10 

3 

4 x lo 15 

713 

114 

20 

1 . 


* These two solutions have decreasing fixed costs (in order gi’ 
+ These three solutions have increasing fixed costs (in order ; 

** Number of times out-of-kilter algorithm applied. 




Augmenting the Algorithm to Consider Additional Side Condit 

There are several additional constraints that might also 
added to the problem to give more generality. Consider the fol 
conditions and the way the algorithm must be adjusted to allow 
their inclusion. 

1. Mutual exclusivity constraints of the form of 

£ y t $ i (2 

Mr E 

where E is some subset of facilities for which only one member < 
subset may be constructed. This constraint is used in the case 
the set E represents several treatment alternatives that might 1 
built at a possible site. Then either none or one of the poasil 
alternatives may be built. It might also be used as an informa 
budget or political constraint. If one facility is built at oni 
site, another cannot be built at a different site. The way tha 
this is handled in the algorithm is as an additional feasibilit; 
requirement. In order to be an exact feasible solution to Prob 
One, the solution to Problem Two must not only satisfy flow thri 
facility constraints but the mutual exclusivity constraints as i 
If it is violated, an exact solution has not been found. A fas* 
convenient way to set up a branching rule for this procedure is 
following: Suppose that the variable chosen for branching i8 i 
and i € E. Since no other member of E may appear in a feasible 
tion if y^* is set to one, set y^* *» 1 and y^ * 0 for all i 4 E 
not equal to i*. Switch all indices i 6 E from S to S. 



2. Budget constraints of the form 


E a^y . £ B (2-4: 

i * W 1 

where W may be a subset of the facilities or the entire set. The 
mutual exclusivity constraint in (2-42) is a special case of (2-43' 
(with all a^ - 1, and B * 1). Each time a member of the set E is i 
candidate for branching, a check is made to see if establishing tin 
facility will violate the budget constraint. If it will not, then 
branching may proceed. If it will, the variable is only allowed li 
the solution at a value y^ “ 0, and the search for another variabl' 
on which to branch is carried out. 

3. Upper and lower bounds on flow through a facility if it 
exists. Equipment may exist for which a certain minimum use is 
necessary to keep it from deteriorating. The constraint takes the 
form 


1 P -r -i 

V n & I *ki « Qi u y 4 < 2 -4' 

k=l 

where and Q^ u are respectively the lower and upper bound on th' 
capacity of the facility if it exists. The change in the algorithj 
that will allow this condition is in step 4. When Z^, the cost of 
the arc, is defined, the capacity of the arc is defined as well. 



Substitute these rules in step 4. 


Z =V 
L i v i 


and the upper and lower bounds are Q^ u and Q 
if i € S and y^=»l 


Z i = oO 



and the upper and lower bounds are zero 
if i fe S and y^=0 


and the upper and lower bounds are Q t u and z 
if i 6 S 


4. Upper and lower bounds on individual arcs such as 



r u 

G ij 


(2 


may be added without any change in the algorithm since the out- 
kilter algorithm allows such restrictions on arcs. 


CHAPTER III, URGE SCALE SINGLE AND MULTICOMMODITY ROUTING IN 

GIVEN NETWORKS 

Introduction 

In this chapter, the problem of routing a flow of commoditi< 
Erom sources through intermediate points to sinks in given networks 
Ls approached. As with the facility location problem of Chapter I' 
:he assumption is that the capacity of the vehicle which comprises 
:he means of transport is smaller than the quantity of material to 
>e moved, so that the small scale collection problems of routing 
jetween sources for collection, or between sinks for distribution, 
nay be avoided. In the narrowest sense, finding the flow through e 
jiven network is a simplification of the facility location problem 
md may be readily solved by the use of a transportation or trans- 
hipment algorithm. However, there are often additional considera- 
tions that might be introduced to the problem which increase the 
-enerality of the formulation and the difficulty of finding solu- 
:ions. The first such consideration is the addition of time and 
:apacity constraints on vehicles, budget constraints on certain 
types of expenditures and constraints on manpower allocation. Work 
.n the general area has been done by Garvin et al. (1957) and Barter 
1967), and work in the specific area of solid waste collection has 


al commodities . The work of Szwarc (1967) on multi- 
anspor tation problems has been extended in this 
and will be discussed. Another area of interest is 
nsidored by Ford (1958), Bomberault and White (1968) 
968) called the dynamic trans-shipment problem, in 
desired to route vehicles not only in a spatial sense 
poral sense as well. This work finds application in 
ng of empty railroad cars and may have application in 
collec tion. 



A. Literature Review 


The earliest report in the literature of a vehic] 
problem through a given network with added conditions * 
Garvin et al. (1957) in an article showing applications 
programming in the oil industry. Included in the probl* 
was that of routing vehicles from bulk refineries to se: 
stations. The problem statement is as follows: 

There is a given set of stations, k, each with a 
product D^, and located at known distances from the bull 
and from all other stations. There are several differer 
trucks available for delivering to the stations, each wi 
C s , a given unit cost of operation and a given number oi 
time period when that type of truck may be used. The pi 
find the delivery schedule that will meet the demands f< 
at the stations within the capacity and time constraints 
minimizing the transportation costs. 

For this particular problem, it is assumed that t 
demands are greater than any of the vehicle capacities, 
presented a linear programming formulation as follows: 


Minimize s 


m 

E 

i=l 


n 

I 

j-l 
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E 

s=*l 


c ijs x i js 
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i=l 
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yijk s 
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u=l 


y^uk 


fc-l, 2, t 

j=l, 2, r 
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subject to: 
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. . . , n 
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k=l 
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J-l 

y ojk " , * D k 
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1 = 1 

x ijs 

= ^ X 

u-l jus 

J-l, 

S®1, 
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(3-5) 
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Z 

k=l 

*ijk 

p 

< Z c s x ijs 
s=l J 

i-1, 

J-l, 

2, . 
2, . 
but 

• * . , n 

. , n 

k i j 

(3-6) 
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n 






z 

i=l 

Z 

J-l 

h ijs x ijs < h s 

s=l, 

2j . 

. . . , p 

(3-7) 




are non-negative integers 


(3-7a) 


c ijs “ cost per unit for shipping from i to j in vehicle type 
x ijs ” the number of vehicles of type s sent from i to j 

« the quantity shipped from i to j with final destination 
= quantity demanded at k 
C g = carrying capacity of s*th type vehicle 
hij s « time required to go from i to j in vehicle type s 

h 8 ■ time that an s type vehicle can be used 

n ® number of customers (or cities) 

p = number of vehicle types 

Expression (3-1) specifies the minimization of the cost of 



vehicle operation. Equations (3-2), (3-3) and (3-4) ar 
trans-shipment constraints and exhibit an incidence mat 
Equation (3-2) specifies that all material shipped into 
destined for sink k must equal the amount of material s 
of the node destined for sink k. Equation (3-3) is a d 
constraint which sets the demand at the k*th station, a 
(3-4) is a source constraint specifying the amount to b 
Equation (3-5) expresses a node constraint on trucks hy 
that the number of trucks of type s that enter a static 
the number that leave it. Inequality (3-6) relates mat 
shipment to the capacity available to ship it, while in 
specifies a time constraint on the availability of diff 
of vehicles. 

The authors have ignored the integer constraints 

solved the problem by using a linear programming code w 

of variables that take non-integer values. The probleir 

becomes unwieldy for any reasonable number of sources a 

For n sinks and s different types of service vehicles t 

2 

2n + ns constraints 

and 

3 2 

n + n s variables . 

A 10 sink problem has more than 200 constraints and 1 1C 
The authors did note the network structure of part of t 
set and suggested that some improvement in running time 
gained by exploiting it. So far, however, nothing furt 



ared in the literature. The linear programming approach con- 
es to appear as bigger, faster codes become available. Barton 
7), in work on aircraft scheduling and routing, wrote essentially 
same formulation and used the same technique to solve it. In 
case, C.E.I.R.'s LP/90 code and an IBM 7094 allowed him to 
e a 202 row, 515 column problem in 16 minutes. The solution is 
rted to have non-integer elements which are rounded to integer. 

In application of models of this form to problems of solid 
e collection, the sources are the locations from which solid 
e is to be collected, the sinks are disposal points and the 
mediate points are transfer stations where some form of pro- 
ing such as incineration or compacting might take place in 
tion to the transfer operation. Arcs joining these points 
esent alternative routes or methods by which the material may 
from sources to sinks. Solving the problem will show the over- 
cost of collection and which alternative routes are used to 
in it. Some authors have looked at the problem of modeling the 
all collection system in this manner. Wolfe and Zinn (1967) 
rted to linear programming to solve the problem for a small 
pie problem. The model they proposed would be an ordinary trans- 
ient model which might then be solved by a primal -dual algorithm 
kly, except for problems that arise in the representation of 
aerators. Any waste flow that enters this node is decomposed 
several different forms whose amounts are given proportions of 
inflow into the node. This means that flows leaving the node 



have coefficients between zero and one and an incidence 


longer exists. Since there are advantages to using a p 
algorithm over a simplex algorithm, Anderson (1968) pro 
of-kilter algorithm for solving this probLem which alio 
nodes to be separated into given proportions. He shows 
examples solved by hand and did not indicate if the alg 
been coded. Thus, while it would be expected that such 
would be much faster than solving the same problem by a 
linear programming code, there is no evidence yet to su 
In many cases in vehicle routing, it is desired 
the use of a common vehicle fleet to move more than one 
between supply and demand points for a commodity. This 
larly true in solid waste collection, since many differ 
of wastes such as residential, industrial, construction 
materials, ashes, leaves and even Christmas trees are p 
collected separately. Adding multicommodity aspects to 
formulations as Garvin et al. (1957) proposed for singl 
flow may be done quite simply. Additional constraints 
shipment network requirements for the additional commod 
the commodities interact in their competition for shipn: 
are easily formulated and added to the problem. Howeve 
formulation is already so complex that even a small prc 
unmanageable in terms of computation technique and time 
of more commodities makes the formulation unworkable. 

To be able to deal efficiently with a tnulticommc 



■m it will be necessary to drop all consideration of additional 
•aints such as those on time, vehicles and work rules which 
•e the simple network structure of the problem. This leads to 
trmulation of multicommodity transportation or trans-shipment 
:ms and to solution methods which will prove faster than linear 
imming using the simplex method. Szwarc (1967) formulated and 
ited a solution technique for a multicommodity transportation 
to which he calls the truck assignment problem. He formulated 
the following manner: 

There exist m production sites each capable of producing p 
ities, and n demand locations where there is a demand for each 

commodities. The demand for a commodity is such that the 

capacity of a vehicle is allocated to only one commodity and 
tined for one demand location only. Thus the problem of 
g a vehicle between supply locations or demand locations is 
d and a vehicle's assignment will consist of loading a corarao- 
t a production source, delivering its contents to one demand 
on and returning empty to a production source for reassignment, 
jective is to find the vehicle assignments that will minimize 
tal distance traveled while meeting the supply-demand constraints. 

This problem, because of the requirements that a vehicle 
ing from source to sink may carry only one commodity and that 
turn from the sink to the source is done empty, is a far easier, 
omplicated problem than the general multicommodity flow problems 
ered by Ford and Fulkerson (1958). Each arc in the traditional 
ommodity flow problem may carry many different commodities at 


the same time. But with the special, formulation, eac 
source of a commodity to a sink of a commodity may ca 
commodity. Thus the problem may be decomposed into a 
simpler single commodity problems in which loaded vet 
shipped from sources to sinks, and empty vehicles frc 
to sources. Arcs representing flow between sources a 
different commodities are not allowed. However, arcs 
back to the sources represent empty trucks and any sc 
reached from any sink. 

Szwarc presents no mathematical formulation oi 
but proceeds to a solution method based on solving tv 
transportation problems for an enlarged network in wi 
of several commodities becomes a distinct source for 
commodities . Sinks are treated in the same manner, 
production sources each supplying p commodities and r 
then he solves first an mp by np transportation prob' 
routing full trucks from sources to demand sites, wit 
tween supply and demand for different commodities blc 
Then, the problem of routing the empty vehicles back 
is presented as another transportation problem. A pi 
shown which reconciles the two solutions into a set < 
assignments. A proof of optimality is included in tl 
is a very quick, efficient method since solution timi 
transportation problems is quite small. A later sec 
chapter gives an extension to multicommodity trans-sl 



the type described by Szwarc with the addition of many 
nsiderations which make the formulation more general. 

Dynamic trans-shipment networks deal with the problem of th 
vement of goods and vehicles from location to location over time 
aling with networks where there are time as well as spatial 
pects to the problem often changes the type of analysis performe 
e concern may shift from routing to finding the maximum flow 
rough a point that may occur in a given time sequence. The 
uting of vehicles must specify not only which arcs of the networ 
e traveled, but in what time period as well. 

The first work to appear in the literature on dynamic flow 
oblems was by Ford (1958) and is also reported in the textbook b 

rd and Fulkerson (1962). The problem was to find the maximum 

namic flow through a given network where the dynamic flow at a 
int is defined as the total flow that passes that point during a 
ecified time period. The normal concept of a capacity of an arc 
a static or time-independent problem becomes a capacity per uni 

me in the dynamic problem and the cost on the arc represents the 

avel time through it. Ford showed that there is a direct relati 
ip between static and dynamic problems and that the maximum flow 
rough the network as obtained by a trans-shipment algorithm coul 
used directly to construct the solution to the dynamic problem, 
mberault and White (1968) were interested in the problem of rout 
pty shipping containers from sources where they exist to demand 
r them. In this case demand exists in only specified time peric 


and there is both a travel time and a cost associated with t 
along an arc. The authors wished to apply a travel cost to 
between nodes as well as a time constraint and then minimize 
cost of shipping the containers to meet the temporal demand 
them. By enumerating the feasible routes that exist for a c 
tainer to reach a demand point before or as it is needed, th 
construct an expanded network through which flows which mini 
transportation costs may be found by the use of a trans-ship 
algorithm, much in the spirit of the Ford and Fulkerson algc 
White (1969) suggested an improved algorithm for solving the 
container problem, which involves decomposing into a nested 
of subproblems and uses an out-of-kilter algorithm for solut 
Applications of the work on dynamic trans-shipment ne 
to problems of solid waste collection are of interest partic 
in the scheduling of large scale vehicles used for transport 
transfer sites to disposal areas. Such vehicles may be eith 
or railroad cars, and questions such as how many of what typ 
needed may be approached using this formulation. 



ormulation and Solution of a General Multi commodity Truck 


ssignment Problem 

Szwarc (1967) described a method for solving multicommodity 
portation problems which is discussed in the preceding section, 
ork is limited in application because it cannot deal with trans- 
ent networks. In the paper on his method, Szwarc did not give 
hematical formulation for the problem he was attempting to 
. However, such a formulation gives a hint as to how it may 
tended to a far more general form, and solved optimally using 
t-of-kilter algorithm. 

Consider the following definitions: 


X. « the number of vehicle loads of commodity k to be sent 
J from production source i to demand site j 

* the number of empty vehicles sent from demand site j 
back to production source i 

c ijk * cost per vehicle load of shipping commodity k 
from i to j 

= cost per vehicle of returning an empty vehicle from 
demand site j to production source i 

Dii. 13 demand in vehicle loads for commodity k at demand 
site j 

S ik * su PP*y in vehicle loads for commodity k at production 
source i 

m “ number of production sources 
n “ number of demand sites 

p « number of commodities 



Then the general mathematical formulation for Szw 


delivery problem is: 
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are non-negative integer 


The objective function in (3-8) is to minimize tb 
loaded trips from source to site and of empty return tri 
Inequality (3-9) requires that the amount of commodity V 
from source i to all sites must not exceed the supply of 
at that source. Similarly, inequality (3-10) specifies 
total amount of commodity k shipped to j from all source 
greater than or equal to the demand for commodity k at t 
Equations (3-11) and (3-12) are flow equations requiring 



: number of empty vehicles sent out of a demand site must equal 
ter of loaded vehicles arriving there, and that the number of 
vehicles sent out of a source must equal the number of empty 
arriving there. Inherent in this formulation is that 
; must return to a source after leaving a demand site. The 
requirement (3-12a) will be shown to be automatically satis- 
ice the problem may be set up in network form and the 
i found by finding the minimum cost-maximum flow using the 
liter algorithm. Figure 3-1 shows the graph corresponding 
s problem for the case of m=2 sources, n=2 sinks, p=2 conuno- 

,e computational difference between the two methods is shown 
‘act that the Szwarc method requires the solution of two HxU 
'tation problems and some manipulation with the solution, 
ilution by the out-of-kilter algorithm a graph of np + mp + 
i arcs and (m+n)(p+l) nodes (or, for the example case, 
and 12 nodes), must be solved, which will take longer to 

There are several advantages of the out-of-kilter approach 
ke it a recommended procedure for general cases, however. 

It is possible to consider any sort of routing network 
ermediate facilities. In this case, it may easily be shown 
ling intermediate node equations to the model formulation 
iaves a graph solvable by the out-of-kilter algorithm, but 
Izwarc’s algorithm. Since most vehicle routing systems, in- 
solid waste systems, involve the consideration of intermediate 
uch an extension is of considerable analytic value. 
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X , x , x , x are non-negative integers 
i j k a j K iak J i 


where : 

a = index relating to intermediate nodes 

i 3 index relating to supply points 

j 3 index relating to demand points 

k » index relating to commodities 

x ° number of truck loads of commodity k sent direc 
J from source i to demand j 



x iak = number of tru ck loads of commodity k sent from 
source i to intermediate point a 

•k 

x ajk 13 nuinber of truc k loads of commodity k sent from 
intermediate point a to demand j 

x,. * number of empty trucks returned from demand j 

directly to source i 



* 

^iak 


= unit cost of supplying demand for commodity k at 

sink j directly from source i = c, + t J( + r 

ijk ik jk 

= unit cost of shipping to a as an intermediate point 
for commodity k from source k = c* ak + fc ik 


^ajk “ unit coSt °f using a as an intermediate point for 

supplying demand for k at j B c ** + u , + r., 

ajk ak jk 

c i1k ” cost P er truck load of commodity k from i 
^ directly to j 

c iak = c ° 8t P er t ruc ^ load of commodity k from i to a 

c a jk - cost per truck load of commodity k from a to j 

Cj^ = cost per empty truck from j to i 

^ik ” co8t P er truc k load of shipping commodity k from 

source i 


r,, * cost per tuuck load of receiving commodity k at 

demand j 

u ak = cost P er truc k load of trans-shipping commodity k 
at intermediate point a 


S ik a supply in truck loads of commodity k at source i 

Dj k « demand in truck loads for commodity k at demand j 

^ijk s u PP er bound on flow of k from i to j 

^iak 53 u PP er bound on flow of k from i to a 

* 

Q a j k * upper bound on flow of k from a to j 

** upper bound on flow of k from j to i 

V Qk » upper bound on trans-shipment of commodity k at 
intermediate point a 



Inequalities (3-16), (3-17) and (3-18) express the constra 
>n (respectively) the amount of each commodity shipped from i, tx 
shipped through a, and received by j. Equations (3-19), (3-20) a 
;3-21) express the conservation of flow constraints that the numt 
if truck loads that leave a point must be equal to the number of 
rucks that enter the point. Inequalities (3-22), (3-23), (3-24) 
,nd (3-25) restrict the flow between points to be less than a boi 
‘he integer requirement (3-25a) is satisfied by the solution of t 
•roblem as an out-of-kilter graph. 

Table 3-1 shows some typical running times for the out-of* 
liter solution using randomly generated data. 


able 3-1. Computation Time for Solving Some Multicommodity Trat 
Shipment Problems Using the Out-of -Kilter Algorithm* 
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0.08 
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59 
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Problem Solver" and an IBM 7094 computer. 
Includes time to generate graph and costs. 


n Algorithm for Constructing a Graph to Represent the Truck 


ssignment Problem 


or each source i, a set of p+1 nodes is created and the nodes 
re named y^, y^, y lp , where y ^ represents the k'th 

ommodity produced at the source i, and y^ represents the 
1 th source. 

irected arcs are drawn from y^ to each of the y^. The lower 
nd upper bounds on capacity of the arcs are zero and S u . 

LK 

he unit cost is t,. . 

ik 

or each intermediate point a, a set of capacitated p nodes 
a = (a^-, a^, a p ) is constructed. The upper and lower 

ounds on flow through the capacitated nodes are V fl j { and zero, 
he unit cost is u a ^. 

or each demand j, a set of p+1 nodes is created and the nodes 
re named qj, q^, qj p , where qj^ represents the k’th 

ommodity received at demand j, and q^ represents the j T th 
emand. Directed arcs are drawn from each node to q^. The 
ower and upper bounds on each arc are zero and D^* The unit 
ost is r . 

irected arcs are drawn from each node y^ to each capacitated 
ode z fl , with upper and lower bounds and zero, and unit 

ost C iak' 

irected arcs are drawn from each node y^ to each node q^> 
ith upper and lower bounds an ^ zero, and unit cost 



7. Directed arcs are drawn from each capacitated node z to e* 


node q., , with upper and lower bounds Q and zero, and ut 
jk r ajk 


cost G 


ajk* 


a. 


Directed arcs are drawn from each node to each node 
upper and lower bounds and zero, and unit cost C^. 



Decoding the Out-of-Kilter Solution of the Graph to O btain 
Routings 

No feasible solution for the graph means no feasible 
solution for the problem. 

The flow from node y\ to node y^ represents the number of 
truck loads of k supplied at i. 

The flow from node y ^ to z q represents the number of truck 

loads of k shipped from i to a. 

The flow from node y^ to q^ represents the number of truck 
loa<^s of k shipped from i to j. 

The flow from node z & to represents the number of truck 

loads of k shipped from a to j. 

The flow from node qj to node y^ represents the number of 
empty trucks returned from j to i. 

The cost of the solution to the graph is the cost of the 
solution to the delivery problem and is optimal. 



collection services. In this chapter, 
vehicle capacity is greater than an it 
it is also greater than the sum of all 
all demand locations before returning 
problems are referred to as single rot. 
route need be fashioned. If the sum c 
collection vehicle must return to its 
completed, and then return to service 
becomes a multi -route problem. 

Referring to demand as discrete 
the context of a network framework. I 
located at the nodes of a network with 
taking place along arcs. Continuous c 
spread along the arcs of the network t 
by having the collection vehicle move 
of a discrete case might be routing pi 
petroleum products to gas stations, wl 
illustrated by snow removal on streeti 
Uniform demand means that the demand : 
the same as any other point such as 
with reasonable accuracy in residen 
uniform demand could be exempli; 
which could require different quanl 

Turning to classifications of 
problems fall into an area of schedul 



CHAPTER IV. VEHICLE SCHEDULING FROM GIVEN LOCATIONS 


Introduction 

The previous chapters have been devoted to large scale "f 
>f goods" problems with the routing of individual vehicles among 
:ol lection tasks neglected. This chapter will deal with the pro 
)f routing of individual tasks through networks for the case whe 
:he capacity of the vehicle is such that it may service more tha 
>ne demand before returning to its terminal. The problems of co 
Lecting from a demand area and delivering to it are analogous an 
Jill be mentioned interchangeably. Such operations are the basi 
service provided by a multitude of agencies in the private and p 
sectors of the economy, and have attracted a great deal of analy 
ittention. 

Distinct areas of application and solution in vehicle 
scheduling may be classified on the following criteria: 

a. What is the capacity of the collection vehicle in 
relation to the demands for service? 

b. Is demand continuous or discrete? 

c. Is the demand for service uniform or nonuniform? 

d. Are there any additional constraints besides requi 
that all demand for service be satisfied by a coll 
tion vehicle that starts from and returns to a ter 

The first main division may be made on the relationship b 
:he capacity of the collection vehicle and the size of the deman 


Lon services. In this chapter, it is assumed that the service 
capacity is greater than an individual service demand. If 
Lso greater than the sum of all demand, the vehicle may visit 
md locations before returning to its terminal. These 
i are referred to as single route problems, since only one 
ied be fashioned. If the sum of demand is such that the 
Lon vehicle must return to its terminal before collection is 
id, and then return to service additional tasks, the problem 
a multi-route problem. 

leferring to demand as discrete or continuous is done within 
:ext of a network framework. Discrete demand implies demand 
at the nodes of a network with travel between the nodes 
>lace along arcs. Continuous demand implies that demand is 
ilong the arcs of the network and is satisfied for any arc 
ig the collection vehicle move along that arc. An example 
icrete case might be routing problems for the delivery of 
tm products to gas stations, while the continuous case is 
ited by snow removal on streets or the delivery of mail, 
demand means that the demand for service at any point is 
: as any other point such as would occur in snowplowing or 
[sonable accuracy in residential solid waste collection. Non- 
demand could be exemplified by the service stations, each of 
>uld require different quantities of product. 

’urning to classifications of problems, discrete, uniform 
i fall into an area of scheduling theory known as traveling 


salesman problems, and research on these problems has been i 
In the traveling salesman problem, a route or schedule must 
up which includes all demand locations once and returns to 
ing point while minimizing the distance traveled. The teri 
traveling salesman problem is generally applied only to the 
route problem, while multi-route problems are called m sale 
traveling salesman problems, where m is the number of route 
fashioned. Since demand is uniform, a knowledge of the cap; 
the salesman or service vehicle and the number of demand po 
be visited, uniquely determines the number of routes. 


I this problem to consider not only that all demands must be 
but that each demand must be serviced within a particular 
'iod or in a certain order. With the addition of such con- 
1 , the single route problem as well as the multi-route 
i may no longer be approached directly as traveling salesman 
Host of the work done for the case without additional 
.nts is adaptable with modification to this more difficult 

If the nature of the demand is uniform and continuous, some 
:ing techniques from graph theory may be brought to bear, 
jle route problems, the Chinese postman problem describes the 
of traveling all arcs of a network at least once in a con- 
tour and minimizing the total distance traveled. The problem 
more difficult if some of the arcs in the network are directed 
0 rather than nondirected (two-way). No evidence of work on 
mte problems in this area has been found, although it would 
:o have great application for many municipal services. 

?urning to the solid waste collection problem, it is obvious 
:ause of the size and complexity of demand it is necessary to 
:h multi-route problems. Arguments may be made as to whether 
tal constraints are applicable. In the next section, a 
ire search of work done in all these areas is presented to 
reeling for the state of the art to this time. Original work 
luthor of this dissertation on some problem areas in m sales- 
>lems is also presented. 


u Literature Review 


1. The Traveling Salesman Problem , This is the name given 
the following problem: A salesman is to visit a group of N ci 
tnd the distances from each city to every other city are known. 
;he tour that the salesman should follow that will allow him to v 
ill the cities and return to his starting point while minimizing 
:otal distance traveled. 

Implicit in the problem formulation is the assumption that 
;hls is a single route problem. This means that the capacity of 
salesman in terms of the number of cities he may visit without re 
:urning to the starting point is greater than, or equal to N. If 
:his is not true, the problem becomes a multi-route problem and i 
ailed the m salesman problem which is defined in the following 
manner: A terminal from which a salesman starts and finishes his 
:ours, and a group of N cities with known intercity distances are 
;iven. Assume that the salesman can only visit k cities (k less 
t) at one time before he must return to the terminal. Find the m 
ours that the salesman must make so that he may visit all the ci 
nd minimize the distance traveled. Each tour must be of exactly 
.ength k plus the terminal, so m is uniquely determined as N divi 
»y k. If the result of the division is non-integer, dummy cities 
lust be added or real cities dropped in order to make m integer. 

By assuming that the salesman must visit exactly k cities on 
iach tour, all considerations of the demand at that city have bee 


nd the possibility of making more than m tours by 
ban k cities on a tour has been ignored, 
irature is full of research on the traveling salesman 
* it contains the elements of so many scheduling 
:h diverse areas as job shop scheduling, aircraft 
.e routing and production planning, but little has 
>out the m salesman problem. Since most solid waste 
>lems are multi-route problems, particular attention 
i this section to those authors* work on the traveling 
im that might be extended to the m salesman problem, 
reling salesman problem, while being rather simple to 
been quite difficult to solve. This is because of 
>r of feasible solutions that exist for even a moderate 
Consider a 20-city traveling salesman problem. For 

L/ 

distance matrix there are 19 factorial possible 
IQ cities may be ordered on the route. This means 
1.2165 k 10 17 possible feasible solutions, which can 
as being a very large number by noting that if one 
computer to enumerate one million of these solutions 
would take over 27,500 years to search all feasible 
is attention has turned towards analytic methods for 
ana to the problem. A definitive and complete review 
Leal development and computational experience with 

atric distance matrix, the distance from city i to city 
sssarily the same as the distance from city j to city i. 


the traveling salesman problem was presented by Bellmore a 
Nemhauser (1968). They included a complete bibliography a 
cuss ion of the literature devoted to the subject. Those i 
in following the development of the traveling salesman prc 
greater detail are directed to that article. In discussir 
approaches to the problem, Bellmore and Nemhauser distingu 
between three basic procedures. These are tour to tour in 
tour building, and subtour elimination. So far, all tour 

y 

improvement methods are heuristics in which rearrangement 
ordering of cities on the tour are made until the analyst 
fied that the ordering obtained is a good solution. At a] 
the tour under consideration is feasible in that each cit) 
only once. Switches in ordering of cities that will imprc 
tour are made, and termination occurs when it appears that 
time has been spent searching or no other improvement sees 
Of course, such methods do not guarantee optimal results, 
do have the feature of being quite fast. 

Tour building involves choosing a starting city anc 
electing which city to visit next. The procedure terminat 
tour returning to the original city and including all othe 

2 / 

- In discussing these methods, the terms heuristic solutic 
optimal solution will be u6ed. A heuristic solution is 
is feasible and that is found by using a set of rules tl 
that it is a "good" solution but there is no assurance 1 
solution found is the best, or optimal. Optimal solutic 
that the procedure, if followed to completion, guarantee 
the answer found is a best solution. 



been generated. An example of a heuristic method using this tecl 
nique is the nearest city rule, which uses as a criterion for 
selecting the next city to be visited, the unvisited city which i 
closest. Tour building has also been used in optimal procedures, 
such as one developed by Little, Murty, Sweeny and Karel (1963). 
particular interest is a technique developed by Pierce and Hatfie 
(1966) to deal with a production sequencing problem which has adc 
nal constraints dealing with job deadlines. The algorithm is bas 
on solving the traveling salesman problem using a tour building 
optimization technique. Reference is made later in this literati 
review to extensions of Pierce's work to single and multi-route 
truck dispatching problems. Held and Karp (1962) developed a 
dynamic programming algorithm for the optimal solution of the 
traveling salesman problem that is limited to about 15 cities by 
storage requirements on present day computers. Some aspects of 
Held and Karp's tour building scheme are also applicable to tour 
building for the m salesman problem but the matter has not been 
pursued In great detail. Such tour building techniques also alw? 
work with feasible solutions. This is not true for subtour eliml 
tion methods, where a less constrained problem is solved which m? 
or may not give feasible solutions. If the solution to the trav? 
salesman problem is feasible, it is also optimal. However, it m? 
be infeasible in the sense that all cities are visited but sever? 
disconnected subtours, instead of one tour, are found. The solul 
is then sought by adding constraints to force the elimination of 


subtours . 


A later section of this dissertation presents a sut 
elimination scheme for approaching optimal solution of th« 
salesman problem, 

2. Single and Multi-Route Truck Dispatching Problems . As 
the preceding section, the problem to be dealt with here 1 
at various times referred to as the truck dispatching pro! 
vehicle scheduling problem, the delivery problem, the clo\ 
problem, the farmer's daughter problem and the truck routj 
Briefly, it deals with the scheduling of vehicles of a gi\ 
among discrete, nonuniform demands for the services afforc 
vehicles. Since demand is nonuniform, the number of rout€ 
fashioned is one of the unknowns. In addition to constrai 
demand for service must be met, additional constraints on 
deliveries, earliest and latest arrival times and vehicle 
are also added by some authors. 

The earliest description of the problem occurs in t 
Garvi'n et al. (1957). In an article on operations researc 
various aspects of the oil refinery industry, the authors 
the large scale problems discussed in Chapter III, and alf 
the problem of routing of individual vehicles. Motivated 
similarity between the vehicle scheduling problem and the 
salesman problem, they named the problem the farmer’s dauj 
problem, since the salesman visited certain nodes more thf 
that presumably was where the daughter was. No attempt at 
tion or solution was reported. 



Dantzig and Ramser (1959) discussed a problem they called t 
ruck dispatching problem. The problem they defined is as follows 
iven N points or cities each with a demand for deliveries of q^, 
nd a terminal point that does not have a demand. Let C = the 
apacity of the vehicle, and in this problem 

N 

max q. < C < £ q . 

i=l 


urther, a distance matrix D“d^ which shows the distance from any 
ity i, or the terminal to any city, is known. The matrix is sym- 
etric (i.e., j=*d j ^ ) . Find the routing for the vehicle that wil 

atisfy the demand for deliveries and the capacity constraint on 
he vehicle while minimizing the distance traveled. 

The authors were unable to provide a model formulation that 
ould solve the problem optimally, but presented an interesting 
euristic which forms the basis for solution of large scale truck 
ispatching problems to this day. This la done by inves tigating t 
roblem where it is desired to aggregate the various demands Into 
airs. This means that for every city i, the vehicle must make 
ither the trip from the terminal to the city or the return trip 
rom the city to the terminal. Since the distance matrix is symme 
he distance traveled by the vehicle between the terminal p, and t 
ities is a constant found by calculating 


N 


E d si' 
1=1 


he problem of minimizing the total distance traveled then reduces 


:o the problem of minimizing only the distance traveled between 
:ities. This may be thought of as finding the pairing of cities 
:hat minimizes the distance between cities of each pair. 


Define x^ * 1 if city i and city j are paired 


= 0 otherwise 


fhen by solving 


N N 


Minimize i z d . .x . . 

i=l j-1 J J 


(4-1 


N 


mbject to: 


L x ij s 1 

i=l J 


j-1, 2, N (4-2 


X ij = *ji 


1=1, 2, ...., N 

j-1, 2, ...., N 


(4-3 


x i j “ (0, 1) 


1 = 1 , 2 , . 

j«l, 2, . 


N 

... N 


(4-2 


'he set of resulting variables gives the desired pairs and cc 


equently the routing for the two city problem. In order to cons 
ime to solve the integer programming problem, the authors sugges 


olving the problem with constraint (4-3a) replaced by 0 ^ x^j «C 


hich may yield fractional solutions in some cases. In that case 
olution is rounded to an integer solution. The problem has beer 
ailed the matching problem, and may be solved optimally very quj 
sing an algorithm by Edmonds (1965). The authors' solution for 
arger problems is to consider the pairs formed at this iteratior 
single point for the next iteration. In this manner routes are 
uilt up until the capacity of the vehicle is met. At each stage 


hors used a simple exchange technique to switch pairings if 
irovement might be found. However, once two cities have been 
at one stage of aggregation, they will never be separated at 
e level of investigation. The authors noted that the method 
istic, and presented computational experience for a 12 city 
\ for which the method does not find what they suspected to 
optimal solution. 

Clarke and Wright (1964) presented a heuristic method based 
:zig and Ramser's technique, but one which will allow vehicles 
.fferent capacities. Their method varies in the means for 
•rating the subparts of the problem and in the means for 
.ng for possible switches at each level. Using the procedure 
iveloped, they found a better feasible solution to the Dantzig 
iser trial problem but could not prove that it is optimal, 
fork serves as the basis for the main computerized vehicle 
Ling package available from IBM as a software service to its 
This package, "System/ 360 Vehicle Scheduling Program", Anon. 

, available for running on the IBM series 360 computers has 
Ln parts, a network analysis section and a schedule producing 
i. The network analysis is designed to find travel distance 
ivel times between potential delivery points and to do some 
Lnary sorting to eliminate obviously dominated schedules, 
iata may be in the form of coordinates or distances and the 
>n of barriers and congested areas may be specified as well, 
gorithm used for finding shortest paths through the network is 


.11s' (1966) work in the decomposition properties of shortest pat 
roblems in large networks. The scheduling procedure has additior 
mristics to allow the inclusion of constraints on such factors £ 
riority ratings, earliest and latest possible starting times, 
iveral commodities carried by one vehicle and different vehicle 
opacities. To date this program is the best large scale analytic 
jol available for analysis of vehicle scheduling problems. 

Balinski and Quandt (1964) presented a formulation for a 
roblem they called the delivery problem, which may be used to sol 
nail multi-route truck dispatching problems. They formulated the 
roblem as a set covering problem which requires that feasible 
Ingle vehicle schedules be enumerated and then chooses from amonj 
lese schedules an optimal set meeting the requirement that all 
smands be satisfied. For a large problem, this would require th« 
great many schedules be generated first before the optimization 
ichnique could be applied. However, if some intuition is availal 
: the problem is highly constrained so that the number of feasib; 
sutes may be reduced, the method offers some appeal. The probler 
tatement is as follows: 

n 

Minimize: 2 c.x 4 (4* 

j-i J j 

n 

2 A.x, « E 

j-i J J 


abject to the constraints: 


(4- 


■j is an m by 1 column vector representing the j'th feasible 
; route, and whose elements 

a^j a 1 if the i'th city is visited on the j’th route 
= 0 otherwise, 

Cj = cost of the j'th feasible route 
Xj B 1 if the j’th route is used 
- 0 otherwise 

E = an m by 1 column vector of l's 
n “ number of feasible routes 
m = number of cities to be visited* 

sasible route is represented as a column vector with a one in 

:h position of the vector if the route visits the i'th city 

;ero if it does not. It is desired to find the set of Xj’s 

.nimize cost and satisfy the requirement that each city be 

i exactly once. Running the problem as a linear programming 

l without integer constraints on the variables does not 

' yield integer solutions. The authors suggested obtaining a 

in using an integer programming cutting plane algorithm. They 

computational experience with a problem of 15 cities and 388 

3/ 

.e routes. The feasible routes were examined for dominance 
tuced to 270. Optimal solution was obtained in 23 iterations, 
ilmilar problem ran 200 iterations without solution. 

•oute dominates another if both deliver to the same cities 
ine has a lower cost than the other. 



The idea of using a two phase algorithm in which the firs 
phase involves generating a set of feasible solutions and the se 
phase is used to select from among that set the optimal routing 
quite intriguing. Some form of branch and bound or selective se 
of the feasible solutions will allow the discarding of many domi 
feasible solutions before the second phase begins. As constrain 
on the charcter of the routes become more demanding, the process 
searching for feasible routes in phase one may even be speeded u 
since many branches may be terminated early and fewer feasible s 
tions become candidates for phase two. What is needed is a bett 
algorithm for solving the set covering problems of phase two, ra 
than cutting plane algorithms, and the work of Pierce provides a 
substantial step in this direction. 

Pierce has written on several aspects of scheduling and t 
dispatching, as well as on the development of combinatorial prog 
ming algorithms for solving set covering problems. As noted ear 
Pierce and Hatfield (1966) made use of a traveling salesman prob 
to look at production problems with additional time constraints. 
This led to Pierce's (1967) major work on truck dispatching. In 
first part of this work the main concern is with single route pr 
with additional constraints that include earliest and latest arr 
times at a specific point, optional deliveries, split deliveries 
constraints on the service vehicles such as volume limit, maximu 
number of stops and maximum time per trip. The solution methods 
suggested by Pierce do not use a set covering approach but are 


ind bound tour building techniques. The basic element on 
ranching takes place is whether or not a particular arc 
from one demand to another is used in the solution, and 
Dund estimates are found using an assignment procedure 
to that of Little et al. (1963). Two ways for searching 
riches are suggested. One is "flooding", in which many 
3 are started at once, and the choice of which branch to 
gate next depends on which seems the most promising at that 
ihe other is to follow one branch to completion {i.e., 
infeasibility, a new feasible solution, or dominance by an 
g feasible solution). Then backtracking may take place up 
ach until it is exhausted and another branch is started, 
ne problem is characterized by extremely long computation 
find and verify optimal solutions. Pierce was more interested 
techniques that generate feasible solutions early so that the 
re may be terminated and the best feasible solution used as a 
ic. For this purpose, he feels that on the average the 
technique will yield feasible solutions sooner. 

In a later paper (Pierce, 1968) the author developed a 
torial algorithm for solving general set covering problems 
se as his trial problems a set of truck dispatching problems, 
tional results based on using the algorithm for different 
oblems are shown in Table 4-1. 

Pierce also found that in dealing with the largest of his 
s, problem 13, the first feasible solution was reached in 


Tab 1 e 4 - 1 . Computational Experience on Truck Dispatching 
of Various Sizes Reported by Pierce (1968) 

Problem Problem Solution Time in Seconds*** for I 

No. Size Forms of the Algorithm 

(m*n)* Alg 1 Alg 2 


1 

5 

X 

31 

0.050 

0.050 

2 

6 

X 

62 

0.167 

0.117 

3 

8 

X 

92 

0.567 

0.200 

4 

13 

X 

91 

35. 167 

6.367 

5 

11 

X 

231 

9.183 

1.383 

6 

11 

X 

561 

12.917 

2.867 

7 

11 

X 

1023 

27.267 

14.383 

8 

11 

X 

1435 

34.950 

19.317 

9 

12 

X 

293 

89.133 

3.500 

10 

12 

X 

538 

131.033 

7.117 

11 

12 

X 

793 

77.767 

4.567 

12 

15 

X 

575 

600.000** 

69.483 

13 

19 

X 

1159 

... 

2400.000' 


m » number of cities, n = number of feasible routes 
Termination without proving optimality 
***\3sing the IBM 7094 




aconds, and the best feasible solution at termination had 
jnd quite early (in 349.767 seconds). In a new work on 
nenta to the combinatorial algorithm, for which only an 
t has been circulated, Pierce and Lasky (1969) indicate that 
;e found an optimal solution to problem 13 in one-half 
This is a considerable time saving, and makes an optimal 
re quite competitive with the heuristic solutions for 

2 size problems for the first time. 

Pierce (1969) is continuing work on solving multi-route 

3 both by tour building branch and bound algorithms, and by 
3e, set covering methods. Garfinkel and Nemhauser (1969) 
eked on a problem of political redistricting using a two 
3et covering approach as well. 

Indrew and Hamann (1967) also suggested a tour generating 
and bound technique. Their work is based on the work of 
at al. (1963), and they branch on whether a particular arc is 
i or excluded from the tour. A single source for starting 
Icles, unequal demands at cities, and route capacities (i.e., 
capacities) that do not have to be equal, are features of 
“1 formulation. A bound on each branch is computed by 
ng the remainder of the possible arcs which have not yet been 
i upon. If a route exists which is below capacity but to 
3 additional load may be added, a branch must be put in back 
starting point. If some branches exist that could possibly 
aute over capacity, they are eliminated as well. Then a 


ower bound may be estimated by using an assignment procedure wh 
akes branch assignments based on minimum cost, but may produce 
ubtours . 

The authors reported a great deal of difficulty in conver 
ence to optimality. They attempted to solve Dantzig and Ramser 
rial problem and found the same solution as Clarke and Wright d 
fter 560 iterations, but could not show that it was optimal. T' 
erminated after 47,000 iterations without having finished the 
otal implicit enumeration. At 14 seconds per thousand iteratio 
ver 10 minutes were spent without having proved optimality on a 
2 city problem. Since a very good solution was developed early 
he authors suggested that the procedure be terminated early and 
he best solution found be used as a heuristic solution. 

Haueman and Gilmour (1967) dealt with a complicated heuri 
ased on solving a traveling salesman problem. They considered 
ingle source and multi-period demands. Their formulation is as 
ollows J 


[inimize 

here 


m 


E C. I Z. + C 2 
j»l i=l J 


D. x j max f . . 7 

( i*Sj ) 
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Dj = min £ 
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(4- 

Z Oj - 1 

j“l, 2, . . . . , n 

(4- 
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E Z ±J - 1 

i“l * 2, m 
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ubject to: 


m 

= w iu ■ 
1=1 

1 

k=l, 2, . 

. . . , m 

(4-10) 

m 

E W . * 

i i ok 

k=l 

n 



(4-11) 

m 

E W io = 
i-1 

n 



(4-12) 

"u = 0 


i-1, 2, . 

. . . , m 

(4-13) 

Z ij ’ M ij 

integers; 

Z ij’ W ik 

>✓ o 

(4-l3a) 


w 

£ ng 1 all a, all j (tour requirements) 

n c s ge j (4-14) 

=» set of all customers in group j 
8 “ any subset of S j ; s“Sj-s 

i « a customer; 0 is the terminal 

» frequency of delivery required by customer i 
d,, a distance from customer i to customer k 

ik 

j « a group of customers, or route; j**l, 2, . n 

n 08 number of groups (a decision variable) 

68 cost per customer delivery 
C 2 « cost per mile of truck travel 
Z i j « 1, if customer i assigned to group j 
« 0, if not 


m 


a number of customers 



w 


ik 


D 


J 


- 1, if truck travels directly from custon 
to customer k 

= 0, if not 

- optimal (minimum) distance to serve groi 


The problem requires in Equation (4-7) that for et 
customers considered the minimum cost route for a vehicle 
This is done using the traveling salesman algorithm deve] 
Held and Karp (1962). For the formulation only the cost 
traveling salesman solution and not the actual routing is 
and the Held and Karp algorithm can produce the cost of t 
quickLy even though finding the route itself takes longei 
Equations (4-8) and (4-9) require that all customers be i 
a group and the terminal be assigned to all groups. Equi 
and (4-11) state that a customer may be visited from onb 
customer, and the vehicle leaving it may go to only one < 
customer. Sending a vehicle back to one self is prohibii 
and the assignment of flow from and to the terminal must 
the number of groups (4-12) and (4-11). Inequality (4-l< 
requirement that states that for any partition of the cuj 
there must be at least one trip of the vehicle between tl 
This requirement is also used in some traveling salesman 
to restrict subtours. For an N city problem there must 1 
constraints . 

To solve the problem the authors generated possib 
of the customers, calculated the value of Dj for each gr 


not directly on the distance between them, but on some indirect 
measure of the compactness of the area. A final collection are 
assignment for a crew per work period does not specify an order 
which it is to be done. The method used for grouping the small 
areas into work assignments is Hess and Weaver's (1965) work on 
political redistricting which is itself a heuristic. The autho 
used this technique to generate several feasible solutions and 
looked for local improvements by switching some assignments. T 
have dealt with a fairly large case of an English city of 150,0 
population and showed a $25,000 per year decrease in collection 
costs for the solution generated over the present method of 
collection. 

3, A Note About the Truck Dispatching Problem . All of the aut 
who have written about the problem to date characterize it as a 
easy problem to formulate and a difficult one to solve optimall 
For this reason, many heuristics have been developed, some of w 
are based on optimization techniques and some of which are not. 
Spitzer (1969) in an article directed towards the practitioners 
scheduling, raised the question of using heuristics to find wor 
solutions versus the use of optimization techniques to find opt 
answers. He pointed out that real world problems may often hav 
advantages over theoretical problems in that they may naturally 
to some form of partitioning or decomposition. While the first 
thing taught in a basic systems course is that the best solutio 


stem is not necessarily the best solution for the system as 
, such approximations might be made along geographic lines, 
the exponential increase in computation time as the number 
=s increase, it is much easier to analyze three 10 city 
3 than one 30 city problem. Hayes (1967) suggested that the 
a of optimization techniques may be to test the heuristics, 
date, little has been done to show how far off a feasible 
a generated by a knowledgeable operator of a complex vehicle 
is from the optimal solution. Pierce (1967) and Andrew and 
(1967) both suggested that the branch and bound routines they 
yised to find optimal solutions are so slow that they should 
as heuristics by stopping the procedure at some specified 
d using the best feasible solution to date, 

Chinese Postman Problem . Turning now to the problems of 
routes for vehicles in networks where demand for service is 
and continuous along every arc in the network, some early 
graph theory paves the way for some interesting analytical 
hes. In 1763, the famous Euler was given the problem of 
a parade through streets of his native city of Konigsberg 
it would cross all seven bridges in the city only once, 
gating the problem, Euler found that because of the configu- 
of the street network and the location of the bridges, the 
s impossible, and in the process postulated some general 
bout finding routes that must cover all arcs of a network, 
r tour is defined as a tour that is continuous, covers every 


arc in a network exactly once and returns to the starting 

4/ 

an undirected network, simple rules exist for checking t 
euler tour exists or not. A node is even if the number o 
touch it is even. An euler tour exists if all of the nod 
If some of the nodes are odd then some arcs will have to 
more than once. In a directed network where all the arcs 
travel limited to only one direction, similarly an euler 
if every node is pseudo-symmetric, or the number of arcs 
towards the node is equal to the number of arcs directed 
it. Rules for establishing an euler tour as well as more 
description of the graph theory involved in the process, 
in the textbook by Berge (1966). 

Finding the optimal route to be taken when an eule 
not exist is called the Chinese postman problem. In this 
route is to be found that minimizes the distance traveled 
more than once. Investigators such as Glover (1967) and 
(1965.) have both suggested optimal algorithms that have l 
successful. Work by Edmonds' associates at the National 
Standards (Santone, 1968), has been proposed for a case c 
waste truck routing in Columbus, Ohio, but the actual woi 
yet been carried out. 



arc. 


a network where travel may be in either directior 



Dptimal Solution Method for the m-Salesman Traveling 
ssman Problem 


in algorithm for the solution of the m salesman traveling 
1 problem under some special conditions is presented in this 
. The problem to be solved is the following: 

Define : 

N = the number of cities to be visited 
excluding the terminals 

s « the number of terminals for salesmen 

k = the number of cities a salesman may visit 
excluding the terminals before returning 
to a terminal 

m « the number of salesmen leaving and entering 
each terminal (the number of tours to be 
established from each terminal) 

5iven, N cities each with the same demand for services. The 
number of cities which a salesman may visit is k. There 
srminals from which salesmen may begin their tours, visit 
k cities and return to a terminal which need not be the one 
Lch the salesman started. Each city is to be visited only 
?ind the route for each salesman that satisfies all of the 
snstraints and minimizes the total distance traveled by all 
ssmen. 

Do eliminate some of the computational problems it is assumed 
;h salesman visits exactly k cities. Thus smk=N, which 
that once s, N and k are known, and N/(k*s) is integer, 


then m is uniquely determined. If N/(k*s) is not an intej 
adjustments such as adding dummy cities or subtracting rei 
are made. 

In the context of the solid waste collection, this 
represents the routing of solid waste collection vehicles 
transfer stations through their individual collection tas' 
traveling salesmen are solid waste collection vehicles, t 
terminals are transfer stations and the cities are small 
areas, each of which generates the same amount of solid w 
time period. The assumption that a vehicle leaves a term 
does not return until it has visited exactly k cities clo 
resembles the actual work rules of collection vehicles wi 
city. Lacking any means of monitoring actual load on the 
during the collection process, the driver is usually give 
assignment and then learns by trial and error at which po 
lection should be stopped in order to dump the material, 
day's work assignment will consist of several collection 
that the assumption that a vehicle does not necessarily r 
the transfer station from which the route started is not 
as long as the last collection route of the day puts the 
at the place from which they started in the morning, in o 
pick up their belongings and transportation home. The as 
that a truck may visit exactly k cities or collection are 
that each collection area contributes 1/k truck loads, I 
this load must be characterized as a load at a discrete p 


it is fairly uniformly spread along the street network to be 
in the collection task. This discrete assumption is less 
} a continuous assumption. As noted in the literature 
continuous problems are Chinese postman problems and no 
ire known for multi-route Chinese postman problems, 

\ important feature of the formulation of the above problem 
Lt will allow more than one terminal from which salesmen 

2 and return. This situation is common to many routing 
including solid waste collection. The work of some authors 

cvering solutions to truck dispatching problems, particularly 
and Quandt ( 1964 ) and Pierce ( 1967 ) might also be readily 
to consider the multi-terminal case. In their work, a two 
ccedure is necessary, with the first phase used for genera- 
3ible routes, and the second phase used for solving the set 
problem necessary to choose the best set of the generated 
routes. Phase one could be adjusted to allow the generation 

3 from multiple terminals and the procedure used to solve 
covering problem in phase two, slightly modified so that 
rminal requirements might be met as well. 


C. A Sketch of the Algorithm for the Solution of the m-S< 


Traveling Salesman Problem 


As described in the preceding section the problem : 
a solution technique is desired is that of finding routes 
exactly k cities for each of m salesmen operating from ea< 
terminals that will visit N cities exactly once and minim: 
total distance traveled, A mathematical programming fortm 
the above problem would indeed be cumbersome with respect 
of variables and constraints necessary to describe even a 
problem. However, the mathematical formulation of a simp; 
which considers all the constraints except that on the nur 
cities each salesman must vi*sit, is capable of fairly stri 
forward exposition. Consider the following: 

Problem B. 


Minimize 


N 

l 

i“l 


N 

T. 

j=l 


ij X ij 


(d 


t=l c=l 


Ct x ct 


a tc x tc) 


subject to * 


N 

l 

c=l 



N 

l 

c=l 



t=l, 2, s 


N 

I 


i = l 


x ij s 


N 

l 

i=l 


i*J 



N 

Z 


i = l 


ij 


j*l 9 2, N 


1 


J“l* 2, N 


N ** 

2 x tc ” m c=1 » 2 > ••••» 8 (4-19) 

c=l 


* -fck 

x ct* x tc’ x ij are non -negative integers (4-19a 


:re 

Xjj = the number of salesmen who travel from city i 
to city j 

x - the number of salesmen who travel from city c 
to terminal t 

x tc = the number of salesmen who travel from 
terminal t to city c 

f i j = distance from city i to city j (d^ =o0) 

d ct = distance from city c to terminal t 

d tc 83 distance from terminal t to city c 

m = number of salesmen dispatched from each terminal 

N = number of cities 

s = number of terminals 

Expression (4-15) is the objective function which is the 
limization of the distance traveled by the salesmen. Equations 
*16) and (4-17) express continuity of flow requirements for each 
'minal and for each city, respectively. They require that the 
nber of salesmen entering a terminal or city must equal the numbe 
iving it. In Equation (4-18) the requirement that exactly one 
Lesman must visit each city is expressed, while in Equation (4-19 
salesmen must visit each terminal. 

No explicit constraint has been written on the number of 


:ities a salesman roust visit before returning to a terminal. Den< 

ls Problem A, the form of Problem B with the added constraint on i 

lumber of cities a salesman must visit. If Problem B is solved ai 
:he solution is a feasible solution for Problem A as well, then ai 
)ptimal solution has been found for Problem A. If not, a branch 
md bound scheme based on solution of Problem B with added 
constraints will be used to find Problem A. 

The constraint set of Problem B has an incidence matrix wii 
i plus one and a minus one in each column, thus indicating it is f 
letwork problem. This insures that solutions will be integer and 

:hat some form of network algorithm may be used to find solutions 

Ln a smaller amount of time than a simplex linear programming 
ipproach. The graph is in circulation form, which means that thei 
Ls continuous flow in the graph from sources to sinks and then ba< 
:hrough directed arcs to the sources. The method of solution will 
>e the out-of -kilter algorithm, but in order to use it some change 
Ln the form of the problem will be made. Equations (4-18) and (4- 
ixpress the requirements that flows through nodes representing 
cities must be exactly one, and through nodes representing termini 
exactly m. This may be accomplished in a graph by the use of a 
capacitated node, as described in Chapter II. 

Each city is represented by a capacitated node with the upf 
cound u^, equal to the lower bound 1^, equal to one. Entering t! 
receptor of the capacitated node are two sets of arcs s set A^ 
represents arcs from other cities to i, and set A* represents arcs 


terminals to i. Note that the empty set. 

he transmitter node of the capacitated node are arcs which 


rise 

two sets : 


is the set of arcs 

which 

go 

to 

other 

om i, 

■k 

and is 

the 

set of arcs that 

return 

to 

the 

termi- 

i. 

B t r\ B* « ( 








a similar manner the terminals may be represented by a 

ed node with upper bound equal to lower bound equal to the 

salesmen who must enter and leave a terminal m. The 

of the capacitated node receives arcs from all cities and 

raitter sends arcs to all cities. 

nsider the following sets of arcs: 

P = (the set of arcs that makes up the capacitated 
nodes representing cities) 

T » (the set of arcs that makes up the capacitated 
nodes representing terminals) 

problem may be shown in simple circulation form ass 
P 

nimize £ 
i“l 


£ X. « 

ifrB (k) 

£ X* k“l, 2, . . ♦ . , a 

(4-20) 

0 * x t 

* 1 i/pUT 

(4-21) 

X i " 

1 i 6 P 

(4-22) 

X i - 

m itT 

(4-23) 

x i - 

non-negative integer 

(4-23a) 


where 


= flow in arc i 
p = number of arcs 

a 3 number of nodes 

C i 3 unit cost of using arc i 


A (k) = 
B (k) = 


0 if i Q P 

0 if 1 fc T 

if i is an arc joining citie 

d ct if i is an arc joining city 
and terminal t 

* 

if i is an arc joining termi 
city c (where these symbols 
defined previously) 

arcs leaving node k 

arcs entering node k 

number of salesmen 


An example of the graph representing this forraulati 
shown in Figure 4-1 for the case of N=4 cities, m=2 salesr 
terminal, s = l terminal and k=2 cities to be visited by eat 
As indicated earlier, such a problem could be solved by a 
algorithm instead of the techniques under discussion now. 
this case is presented here to give an example for a grapl 
easily drawn and understood. The number of arcs and nodej 
to characterize a problem is quite large. For an N city, 
problem when all cities may be reached from all other citj 
terminals, the number of nodes is 


2s + 2N 




(UPPER BOUND, LOWER BOUND, UNIT COST 
NOT ALL POSSIBLE ARCS ARE SHOWN, 


Figure 4-1. Graph for a Problem with N=4, p=l, m=2, k=2. 



and the number of arcs is 


N 2 + 2sxN + s 

For N=16, s=2 there are 34 nodes and 322 arcs. 

The relationship of the solution of Problem B to Pro 
may be summarized in the following three cases: 

Case 1 . The solution of Problem B is also a feasibl 
solution for Problem A. This is shown in Figure 4-2 
This case will be referred to as a feasible solution 
Case 2 . The solution of Problem B is such that m 
routes are found from each of the terminals, but not 
all visit exactly k cities before returning. All 
cities are visited by a route starting from and re- 
turning to a terminal. Such a solution is shown in 
Figure 4-3. This will be called a regular infeasibl 
solution. 

Case 3 . The solution of Problem B generates more th 
ms routes and only ms of these start and end at a te 
A feasible solution, to Problem B might contain some 
that start at a city, visit other cities and then re 
to the starting city without ever touching the termi 
This will be called a phantom infeasible solution, 
example is shown in Figure 4-4. 

If an infeasible solution such as Case 2 or Case 3 a 
additional work will have to be done to find the optimal so 




PPER BOUND, LOWER BOUND, UNIT COST ) 


Figure 4-2, Feasible Solution to Problem A. 






Co Problem A. The type of infeasibility will suggest the way that 
Che branch and bound procedure can be developed. 

The general plan of attack will be the same as for the 
capacitated trans-shipment facility location algorithm. Solve 
Problem B. If the solution to Problem B is Case 1, it is feasible 
to Problem A. Then since Problem A is more highly constrained thi 
B the optimal solution to A has been found. If the solution to B 
of Case 2 or Case 3 and therefore infeasible to A, the problem wi' 
be broken into a series of subproblems which divide the set of 
feasible solutions to B into mutually exclusive collectively 
exhaustive subsets. This is done by taking some infeasible route 
the problem solution and prohibiting arcs involved in that route 
from being in a further solution. Suppose the solution to Problei 
contains a route which starts at terminal t, visits h cities, 
a X> a 2* ••••» a h and retur ^ s to terminal t . Then the set of 
feasible solutions may be divided into h + 1 subsets with the firi 
subproblem having all arcs broken from t to all the cities 
through a^. The second subproblem breaks all arcs from to all 
other cities (a^ through a^) and the terminal t « This is contin 
until the ta. + l'st subproblem breaks all arcs from a^ to throu 

•ft 

a h-l and t0 t * The arcs are broken by setting their unit costs 
infinity. Each subproblem is solved to find its lower bound. Si 
in each subproblem an arc has been broken which previously had fl 
into it, the flow pattern will change and the cost of the new 
solution will be the same or higher than the previous solution. 


cost of each of the subproblems is examined and the one with the 
least cost is found. If its solution is feasible to A, an optirr 
solution to A has been found. If it is not feasible to A, more 
subproblems are generated in the same manner and the lowest cost 
subproblem is again investigated. This is continued until the 
lowest cost subproblem is a feasible solution to Problem A. 

The rules used to set up the subproblems are called break 
rules. Since there are two different types of infeasible routes 
it is possible to have several different break rules, One break 
rule might be to choose for branching the smallest regular 
infeasible route. Another is to choose for branching the smalle 
phantom route and break its arcs. Experience with both break 
rules is reported. 



D. A Detailed Algorithm for the m-Salesman Traveling .Sales 


Problem 

1. The following data are given: 

N = the number of citie3 to be visited excludin 
the terminals 

s = the number of terminals 

m *» the number of salesmen leaving each termina 

k = the number of cities excluding the terminal 
that a salesman must visit before returning 
to a terminal 

The number of salesmen, m, is an integer determined 

If this calculation is not integer, N is adjusted b} 

dummy cities or subtracting real cities until m is : 

d distance from city c to terminal t 
ct 

d. = distance from terminal t to city c 

U V 

= distance from city i to city j 

These distances are not necessarily symmetric. Tha 

* 

need not be equal to d^ a nor f^ equal to f j i • 

2. Construct a graph in the following manner: 

There will be two sets of nodes, P and T, which 

the graph. 

Define P *= (the set of nodes that represent cit 
For each of the N cities, a capacitated node is 
lished. This is accomplished by representing the c 
node by two nodes with a directed arc between them, 


in Figure 4-2. The node that sends the directed arc is 
know as the receptor node, and the node that receives the 
directed arc will be known as the transmitting node of the 
capacitated node. It will be required that the flow througl 
the capacitated node be exactly one unit. This is accom- 
plished by setting the upper and lower bound on flow througi 
the arc joining the receptor and transmitting nodes equal t< 
one unit. Refer to the capacitated node representing a cit; 
as a city node. 

Define T = (the set of nodes that represent terminals) 

For each of the s terminals, a capacitated node is con- 
structed in the same manner as for the cities. It will be 
required that the flow through each capacitated node repre- 
senting each terminal be exactly m units, since m salesmen 
will pass through each terminal. This is accomplished by 
setting the upper and lower bounds on the flow through the 
directed arc joining the receptor and transmitting node as 
m units. Refer to the capacitated node representing each 
terminal as a terminal node. 

Arcs of the graph are defined in the following manner: 

From T to P i define a directed arc from the transmit tii 
node of each terminal node to the receptor node of each cit; 

4U 

node. The unit cost on the arc is d fcc and there is an uppe 
bound of one on the arc. 


From P to T : define a directed arc from the tram 


node of each city node to the receptor node of each t( 
node with a unit cost of d cC and an upper bound of ont 
Within P; from any city node i to any other city 
i ^ define a directed arc with unit cost f^j and at 
bound of one, from the transmitting node of i to the 
receptor node of j. 

3. Define S = (set of all arcs of the graph that ha^ 

excluded from the solution and hence 
have zero flow) 

S = (set of all arcs of the graph not yet 
from the solution) 

S = (set of all arcs of the graph) 

•ff I 

Then Sy S_ = S and S f\ S = , the erapt; 

Initialize S * 

* 

S - S 

4. Find the minimum cost, maximum flow through the g: 
using the out-of-kilter algorithm. If this is the fi: 
step 4 has been entered and there is no feasible solu 
to the graph, there is no feasible solution to the pri 
Define the cost of this flow as CGRAPH. 

5. Examine the solution of the graph to see if it is 
feasible solution to the m salesman traveling salesman 
This is done by decoding the solution to determine th 
for each salesman. From each terminal node only m ar> 


leaving it may have non-zero flow and for each city only one 
arc entering and one arc leaving it may have non-zero flow. 
The routes for each of the salesmen may then be calculated 
by starting at a terminal and selecting an arc to a city, say 
city i , from that terminal with non-zero flow. Then examine 
the city node i , and find the single non-zero arc leaving it 
If this arc goes to a terminal the route is complete. If it 
goes to another city, say i , examine the city node i and 
apply the same rules until the route is terminated, Continue 
until each of the m routes at each of the s terminals has 
been enumerated. If each of the routes enumerated visits 
exactly k cities, the solution is feasible, 

6. If the solution is feasible, go to step 7. If it is not, 
one of the infeasible routes is chosen for branching. Select 
the smallest route in terms of number of cities visited as 
the candidate for branching. In case of a tie for the 
smallest, choose the first enumerated, 

7. Store information about the solution. File on a list, 
ranked on the cost CGRAPH, of all solutions generated, data 
about the feasibility, members of the S set and the candidate 
for branching. 

8. Remove the first solution on the list. If that solution 
is feasible, it is optimal. 



9. If the solution is infeasible, several subproblems 
set up, Suppose the route chosen as the candidate for 
branching starts at terminal t, visits h cities, a^, b 
a^> and returns to terminal t . Then h + 1 subpi 
will be generated. First eliminate all arcs listed ir 
S set by setting their costs to infinity. 

9a. Subproblem 1 - Break all arcs from t to each 
the cities a^, a 2 » . a^, by setting the costs 
these arcs to infinity and add these arcs to the £ 
Repeat steps 4 through 7. Restore the arcs just \ 
to their original costs and go to 9b. 

9b. Subproblem 2 through h + 1 - If this is the { 


subproblem, break all arcs from a to each of 


other h-1 cities and to t and add these arcs to 


S set. Repeat steps 4 through 7 and restore the i 
just broken and go to the next subproblem. When i 
h + 1 subproblems are completed, go to step 8. 


tational Experience 


e algorithm was programmed in ASA Fortran IV and run on an 
computer. The out-of -ki Iter algorithm was written by 
d made available through the IBM SHARE Library. The list 
g subroutines used for keeping track of solutions were 
y Bellmore (1966). Two possible break rules were considered 
ng up subproblems. Break rule A selected the smallest 
e route that passed through a terminal as the route for 
. Break rule B involved selecting a phantom route (Case 3 
lity) which by definition does not pass through a terminals 
ndidate for branching. If no phantom routes exist, rule A 
ed. Both rules were tried in early runs for some 8 city 
and showed little difference. Only break rule A was used 
rger city cases. Because of dimension limitations and the 
cimal packing instead of binary packing for the storage of 
on about solutions, many problems were limited by exceeding 
ather than by time limits. Future attempts at working with 
rithm will require better programming techniques and more 
acking of data to insure greater efficiency, but such pro- 
efforts were not deemed necessary for this trial demon- 
of the method. Runs were made with one and two terminals 
of 8, 12 and 16 cities. Results of these runs are shown 
4-1. Runs were made with symnetric and asymmetric data, 
trie case generally took longer to solve. This same 


property in traveling salesman problems is explained by Bel 
and Malone <1968). The procedure has a tendency to produce 
few feasible solutions before an optimum is found. This wo' 
a considerable disadvantage in the last two runs, since con; 
time was invested before termination without finding one fei 
solution. The algorithm could be speeded up and have more 
solutions generated, by taking each infeasible solution and 
some heuristic procedure to round it to a feasible solution 
development of the feasible solution would not take long, a 
found to be lower than the current best feasible solution, 
bound on the problem has been found. At termination before 
solution, a good feasible solution would still exist. 

Edmonds and Johnson (1969) are presently working on 
area known as degree constrained subgraphs which would appe 
have application to the above problem. Basically, the prob 
to minimize the cost of using arcs subject to constraints t 
specified number of arcs must enter and leave each node. T 
form the basis for generating tours for each vehicle. The 
Edmonds and Johnson deal with are as fast as the matching p 
algorithms. Based on Bellmore and Malone's (1968) experien 
would be expected that this would alleviate the problem of 
different computation time for symmetric and nonsymraetric d 


cases. 


-Cf to 

M 0 O' d) 


I \J 'I-' r-* V* 


o -a o « 

5a eh o 


<0 <u 


4J p- 

a# 


4h w 

o <u 


u 

•H '■*> 


. +J JS 


O -H ^ 

, 

5a o 



t 


<t 

V) rH 

4> 

4h 0) 43 tO 

rH 

O rH U P 

X 

0 Cd *3 

BJ 

. -H <U B 


0 X H 


z 0) U Q) 


cn 

o 

CO 

m 

PH 

CN 

CN 

rs 

On 

* ■ 
r- 

Jc - 

CM 

k 

rH 

o 

r-H 

<1- 

f-v 

m 

m 

cn 

00 

cn 

CO 

rH 

CN 

o 

o 

o 

o 

o 

o 

O 

T— * 

CN 

rH 

cn 

cn 


CN CN ri N N N 


<r co cn o o 


o 

rH 

rH 

NO 

rl 

o 

in 

O 

O 

to 

in 

o 

o 

Nt 

in 

rH 

rH 

cn 

n 

<u 

rH 

rH 

rH 


CN 

CN 

CN 

CN 

cn 


o t" 
co cN 
X 


o 


co 


o 

o 

o 

o 

o 

o 

o 

o 

lO 

00 

CO 

Nf 

CN 

co 

o 

CN 

o 

o 

CN 

00 

o> 

m 

\D 

o 

f- 

o 

1 — l 

*. 

* 

A 

#■ 

A 

o 

•> 

* 

X 

un 

rH 

LA 

r-» 

cn 

o 

Nt 

m 

Nf 

ID 

iO 

tH 


NT 

<J- 

m 

rH 

CN 



coco VOVO'0'0'D-J'D^^4 

^-IrH NC'lNNNf’llOOtO'J 


sRSSSSssssSs 


rH rH 


1 1 1 II s B 11 1 g| 




CN 


0] CN <t <T 


N FI tO N 4 


<J- CN •sJ' ^ 


> « H 



r-l *“"* 


rH H N 


N N N N 


X 

*rt 

(0 

« 

4) 

4h 


T3 

<u 

•w 

« 

u 

0 ) 

G 

<u 

00 


Tt 
G 
3 
O 
td 


rH 0 
O 01 


<U 


3 X 
O 'H 
X M 
td 

U 4> 

d) 4H 

8 9 

H G 


X 

X 


jJ 

p 


•<H 

*r-l 



3 


0» 

4) 

cd 

t>0 

00 

44 

td 

cd 

<d 

u 

a 

"O 

0 

o 


4J 

*j 

O 

(0 

w 

•H 


p 

44 

44 

X 

O 

O 


a) 

a) 

1 

w 

« 


G 

G 

« 

cd 

<d 

« 

0 

V 


4) 

<0 

w 

X 

X 

G 

td 



CHAPTER V. THE ANALYSIS OF A SOLID WASTE COLLECTION SYSTEM 


Introduction 

In this chapter, a large scale solid waste collection am 
disposal system will be investigated, using some of the techniqi 
suggested in previous chapters. The system chosen for examinat: 
was that operated by the City of Baltimore, Maryland. This 
selection was based on several factors. First, Baltimore is a 
good example of a large city with extensive investment in a pub 
solid waste collection system. Second, a great deal of data hai 
collected concerning the operation of the system, and is avail. 
Analysis of the Baltimore system has already been attempted by 
Truitt (1968, 1969), who used the city as an example in the 
development of a simulation model of collection. Through his 
excellent work, considerable data and insight into the system a 
available for use with the models proposed in this study. All 
taken from Truitt's work will be referenced. 


acfcerization of the Baltimore Solid Waste Collection System 


he City of Baltimore, Maryland, had a population in 1960 of 
r one million, a land area of 80,3 square miles and physical 
r is tics that ranged from high density urban slums to verdant 
ity residential subdivisions. The municipal functions of 
ste collection and disposal as well as general sanitation 
es are carried out by the Bureau of Sanitation of the 
nt of Public Works. For purposes of logistics and super- 
the city is subdivided into five autonomous districts, the 
tern, the northeastern, the central, the eastern and the 

These districts are shown on the map in Figure 5-1. The 

1 / 

mploya 1271- persons and owns approximately 340 vehicles 
es of equipment, and operates on a budget of over seven 
dollars per year. A breakdown of the functions and their 
shown in Table 5-1. As indicated in the table, collection 
for almost 85 percent of the department's operation and 
ominant collection operations are mixed refuse collection 
et cleaning. 

he Bureau collects mixed refuse from residential and non- 
al sites twice a week, once on either Monday, Tuesday or 
y, and once on either Thursday, Friday or Saturday. Since 
on in the earlier part of the week represents the pick up 
days 1 accumulation of wastes, crew size is one driver and 

eport of Department of Public Works (Anon, 1966a) 



three laborers. In the later part of the week crew size is re 
to one driver and two laborers. Two types of vehicles are gen 
used for collection. They are both equipped with compaction e 
raent and have capacities of 20 and 13 cubic yards of compacted 
wastes each. The larger vehicle is often referred to as a fiv 
truck and the smaller as a three ton truck. The larger vehicl 
thought to be more efficient than the smaller and is assigned 
general tasks while the smaller is saved for alleys and other 
restricting jobs. In 1965 the Bureau collected 336,893 tons o 
mixed refuse, 31,395 cubic yards of ashes and 216,622 cubic ya 
street dirt. The collection of the mixed refuse from residenc 
involves 92 vehicles, each with three daily route assignments 
week for a total of 276 routes. There are two city owned inci 
tors for the disposal of wastes, and the solid waste collected 
taken directly from the collection area to the incinerator by 
collection vehicle. One of the incinerators is the Number Thr 
(Reedbird) with a capacity of 600 tons per 24 hour day, and th 
other is the Number Four (Pulaski) with a capacity of 800 tons 
24 hour day. The location of these facilities is also shown i 
Figure 5-1. 


1. Cost of Operations and Services Provided by the Bureau 
of Sanitation, City of Baltimore, Maryland, 1965. 


Operation 

Yearly Expense 

Percent of Total 

et Cleaning 

2,792,000 

37.2 

Cleaning 

202,000 

2.7 

d Refuse Collection 

3,034,000 

40.4 

r Miscellaneous 

324,000 

4.4 

ion Expenses 

6,352,000 

84.7 

aerators 

1,110,000 

14.8 

s 

40,000 

0.5 

1 Expenses 

1,150,000 

15.3 


7,502,000 


100.0 


B. Study Objectives 


Analysts involved in the study of, and planning 
scale public systems such as the solid waste colleetioi 
systems must be concerned with both long range and shoi 
planning. Not only must they study the evolution of m 
and methods that may dramatically change alternatives > 
in the future, but they must be concerned with the prei 
vital existence of the system. This study will be dir< 
the analysis of the system within the short range pictt 
to give some feeling for the alternatives that are ava 
immediate introduction into the system. These alterna 
physical changes such as the building of new structure, 
changes such as changes in the manner and character of 
and rules of operation of the system, or changes in th 
carrying out certain objectives of the system. Short 
tives that fall into each of these categories make up 
issues that require analytical attention. 

The three questions that make up the theme of t 

1. Are transfer sites within the city a feasib 
If so, where should they be located and to 
should they be built? 

2. What is the cost and effect of increasing c 
frequency from twice a week to three times 

3. Under what conditions and at what price wou 
haul become a feasible alternative for the 


Subordinate to these main questions but still of great 


lyst and planner are such additional questions as: 


4. How sensitive is the system to changes in parameters 
and cost estimates? 

5. What are the effects of political and esthetic 
constraints that might force a change from the 
solution most economically efficient from the 
engineering point of view to one that is perhaps 
more acceptable to segments of the community? 

6. Are there advantages to cooperation between 
governmental units within a region where such 
cooperation does not now exist? 


uestions will be subject to some preliminary investigation 
next sections, using the large scale facility location model 
ter II. More detailed study would also require some of the 
techniques of Chapter IV to achieve overall efficiency of 
on, but these have not been included in this overview. 


C. Calculations and Data for the Baltimore Stud~v 


■ 


l t Collection area to be Investigated 

The area chosen for study was the same area sti 
Truitt (1968, 1969). It comprises the northwest divis 
City of Baltimore, and a map of the area and its relatj 
the city are shown in Figure 5-1. Within this area it 
to subdivide into smaller subareas for analysis. A cc 
to do this is to consider each census tract within th£ 
subarea. Exact population and neighborhood density t) 
the 1960 census , are accurately known. There is a tot 
census tracts within the northwest quadrant and these 
the map of Figure 5-3, A coordinate system has been £ 
zero point at the northwest corner. Information about 
number of housing units, neighborhood density type am 
coordinates of the approximate centroid of the tract t 
Table 5-2. Neighborhood density types are based on et 
from census data by the City of Baltimore, Department 
as reported by Truitt (1968). These definitions are £ 
Table 5-3. The housing unit density is important bec£ 
some indication of the speed of collection within a si 
would be suspected that the more closely spaced the pj 
are, the faster the collection process may be in ternu 
per hour collected. Truitt, in a statistical study oi 
ata, shows that there is a significant difference bel 


able 5-2. Pertinent Data Concerning the Forty Census Tracts of th 
Northwestern Division, Baltimore, Maryland 


umber 


Neighborhood 

Density 


Number of 
Househo ld 
Units 


Population 

1960 


Coordinates 



Tract number is U. S. Census designation. census. 

Population and number of household units fr0 * ce ntroid o 

Coordinates are in feet measured to the appr nor thvest cor 

tract. Point x-0, y=0 is located at the extreme nor 

of the Baltimore City boundary. 


type 1 and density types 2, 3, and 4, but no signific 
between 2, 3, and 4. Thus only two sets of neighbor! 
types need be considered, type 1 and other than type 


Table 5-3. Classification of Neighborhood 

Classification Housing Units per 

1 Ten or less 

2 
3 


4 


Eleven to twenty 
Twenty -one to for 
More than forty 



5-1. Operating Divisions of Bureau of Sanitation, 
Baltimore, Maryland. 





Figure 5-2 


Map Showing Census Tracts in Northwe 
of Baltimore, Maryland, 


cgro 



Figure 5-3, Hap Showing Location of Proposed Transfer Sites and 
Present Incinerator Sites, Baltimore, Maryland. 


2. Speed of collection within a subarea 


The speed of collection within a subarea is a f 
neighborhood density and the number of days since the 
tion. When the collection frequency is two times a we 
Baltimore, one collection is of four days of solid was 
and the other of three. Normally the City of Baltimor 
driver pluB three laborers for the larger collection, 
plus two laborers for the smaller one. 

On pages 174-176 of Truitt (1968) several histo 
presented for collection rates in pounds per hour as a 
days since last collection and neighborhood type. The 
based on twice a week collection are computed from act 
tions within the system. Estimates were made of histo 
times a week collection frequency. In this case the f 
of the week has three days accumulation and the remain 
two days accumulation. The number of personnel assume 
and two laborers. Table 5-4 presents weighted average 
grams for different conditions which will be used in t 
models . 

3. Cost of collection within a subarea 

Cost of collection is assumed to be a function 
frequency, the speed of collection, the neighborhood t 
tion, and the work rules assumed. The following data 
to this computation: 



i-4. Average Collection Rate in Pounds per Hour for Different 
Days Since Last Collection and Neighborhood Type 


nee Last Neighborhood Type Average Pounds per 

ction Hour Collected 


2 1 2,500 

2 2, 3, 4 4,580 

3 1 3,000 

3 2, 3, 4 6,300 

4 1 5,000 

4 2, 3, 4 6,300 


Table 5 - 5 . Collection Costs in Dollars per Hour 


nee Labor Necessary 

ion 

Labor Cost 

Do liars /Hr 

Vehicle Cost 
Dollars /Hr 

Total 

Cost/Hr 

1 driver + 

3 laborers 

12.63 

4.40 

16.73 

1 driver 4- 

2 laborers 

9.33 

4.40 

13.73 

1 driver + 

2 laborers 

9.33 

4.40 

13.73 


a) Collection vehicle costs - The City of Baltin 
estimates its cost per hour of vehicle use for the 20 ct 
collection vehicle as $4.40 per hour. 

b) Typical labor rates are twenty dollars per ds 
driver and eighteen dollars per day for a laborer. The 
assignment for a collection of four days since last coll 
one driver and three laborers per vehicle. This drops t 

and two laborers for shorter intervals between collect! 

c) To compute an hourly labor rate from the dail 
is necessary to compute the average time spent in actual 
Truitt (1968) has estimated that between 45 and 50 minut 
are spent in actual productive time, with the remaining 
spent in start up, repairs, off service time, etc. For 
of this study, 45 minutes per hour or 6 hours per day wi 
as the productive time. This yields a productive hourly 
drivers of $3.33 per hour, and for laborers of $3.00 per 

Table 5-5 shows labor, vehicle and total costs fc 
types of collection. 

d) Waste load generated - The waste load generat 
person per day is a statistic that shows great variation 
(1968) after extensive study of Baltimore residential da 
1.95 pounds per person per day and calculated a standard 
of 0.09 pounds per day. This figure will be varied to s 
sensitivity of the solution to it. 

e) Total collection cost in a subarea per week - 



llection cost per week in a subarea is calculated in the following 


aner: 


Define 


Pi * population of the i‘th subarea 

k^ = weight of solid waste generated per person 
per day with twice a week collection 
frequency 

k i - weight of solid waste generated per person 
per day with three times a week collection 
frequency 

lj = labor cost per hour for a crew of one 
driver plus j laborers 

®ik = c °H ect i° n rate in pounds per hour for 
neighborhood type k for i days since 
last collection 


The reasons that the waste generation is shown to be different 
different types of collection frequency are pointed out in stafcis 
al studies by Quon, Tanaka and Charnes (1968). They show a 
nificant increase in waste per household per week when collection 
quency increases from once to twice a week. 

Total collection cost in the subarea is C^. 


C i 




l 3 

8 ^» j 



dollars per week for a 
collection frequency of 
twice a week, and 



(3 _i*_ + 2-ii- + 2 -la_) 

g 3,j g 2,j S 2>J 


dollars per week for a 
collection frequency of 
three times per week 


C i 

collection cost per ton = x 2000. 

7 Pi ki 


extra charges were made for overhead or supervision as there was 
way to obtain a reliable estimate of these figures. 



4. Cost of Transport by Collection Vehicle from t 
to Transfer Facilities or Disposal Points 


a) The speed in traffic of a collection vehicl 
not to differ significantly between a loaded and an un 
condition by Truitt (1968), whose histograms also indi 
average speed of 16 miles per hour. 

b) The cost of the vehicle and crew per hour i 
calculated in 3 (c) above for collection within a sub 

c) The distance between the subarea and a disp 
fer site is calculated in the following manner. A set 
coordinates are given for all transfer stations and di 
The location of the subarea is taken as its centroid, 
vehicles usually cannot travel in straight lines betwe 
travel distance is estimated using a metropolitan or r 
measure in which the travel distance between a and b i 



where j | means absolute value. 

d) The average weight carried per 20 cubic yar 
vehicle at the disposal point is indicated in a histog 
in Truitt (1968). This is taken as 9,000 pounds per v 
the cost per truck for travel from the disposal point 
collection area and return is measured in the followir 




Jefine 


d 


- travel distance in miles in one direction 


1^ = cost of vehicle + driver + i laborers per hour 
s = travel speed in miles per hour 
hen cost per truck load for the two way trip is 


2 d 
s 


X 


k 


e each truck carries on the average 9,000 pounds or 4.5 tons, 
per ton is 


_1 

2.25 



1 


i 


a 5-6. X, Y Coordinates for Proposed Transfer Stations* 


nation Type 

>1 

X Coordinate 
(feet) 

Y Coordinate 
(feet) 

Transfer to 
railroad cars 

24,000 

21,000 

Transfer to 
trailer trucks 

16,000 

11,000 

Transfer to 
trai ler trucks 

14,000 

22,000 

Transfer to 
trailer trucks 

28,000 

20,500 

Transfer to 
trailer trucks 

22,500 

7,000 

Transfer to 
trailer trucks 

37 , 000 

9,500 

Transfer to 
trailer trucks 

34,500 

21,000 


, figure 5-3. 



Transfer Facility Costs 


Truitt (1968) listed four potential transfer si 
located within the northwest quadrant, which will be u 
study as well. In addition, three more sites were cho 
quadrant at various locations including one that is in 
transfer to rail-haul selection. The coordinates X, Y 
site are given in Table 5-6 and are shown on the map i 
The original Truitt sites are marked A, B, C, and D. 
additional vehicle transfer sites are E and F, and the 
is RR. The cost of building transfer facilities was a 
by Truitt (1968) in the following manner: 

Land Cost: Fixed cost of $40,000 below 200 

capacity and a variable cost of 
of additional capacity above 200 

Labor Costs: Three men at $20 per day for a t 
station of less than 200 tons pe 
per day. Four men at $20 per da 
day for a greater capacity. 

Utility Cost? $1200 per year. 

Capital Cost 

of Structures: $125,000 for a capacity of less 
per day. Above 100 tons per da^ 
variable cost of $500 per ton of 
capacity added to the fixed chai 



The to Cal waste load per week from the northwest quadran 
estimated as 14 pounds per week for each of the 203,500 persons 
1428 tons per week. The peak day load would occur in the early 
of the week with three times a week collection. In this case o 
Monday and again on Tuesday, one half of the northwest quadrant 
visited and 3/7 of the weekly waste load is collected. Thus th 
maximum daily load is 

1/2 x 3/7 x 1425 = 306 tons per day 

The cost of a transfer facility as a function of weekly 
capacity is shown in Table 5-7. The yearly cost of the structui 
has been produced through discounting, assuming 30 year life anc 
interest rate of 10 percent. The yearly cost of the land is bas 
on interest on investment and it is assumed that the land does r 
depreciate. No estimates were made on taxes lost on the land, 
since Truitt's sites A, B, C and D were chosen on land that was 
already public. It is assumed for simplicity that the same type 
of land is available at E and F. 

6. Transfer from Transfer Station to Final Disposal Via 

Transport Vehicles and Disposal Costs 

a) Transfer costs - The vehicles suggested for transport 
the final disposal areas from the transfer station are of two ty 
One type is a 75 cubic yard tractor trailer with a weight capaci 
of an estimated 35,000 pounds. Truitt (1968) gives the cost of 


Table 5-7. Cost of Transfer Stations as a Function of Weekly Capacity 
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le at $11.00 per hour. Using $3.33 per hour for a driver gives 
t of $14,33 per hour of truck use. The average speed is esti- 
at 16 miles per hour. The other type would be rail-haul from 
ransfer station to some far distant point for disposal. 
ielphia, Pennsylvania, as reported by PUBLIC WORKS magazine, 

, 1968), presently does this at a contracted rate of $5.40 
m for shipping and disposal in abandoned coal mines 200 miles 
.n western Pennsylvania. This cost will be used as an 
ite of rail-haul transport and disposal costs. 

b) Estimates of tractor trailer vehicle cost per ton. 

Define Xp, Y^ = coordinates of the facility 

X^, Yq = coordinates of the disposal point 
c 53 capacity of the vehicle in pounds 
s = travel speed in miles per hour 
d = distance to be traveled - one way 
1 - cost per hour of vehicle plus labor 

ost per ton => (2 - x I x 2000) 
s c 

d is measured as before by a rectangular metric 

d = fx F - x D | + |y f - Y D | 

| | means absolute value of the difference. 

c) Costs and locations of disposal areas - The location of 
sposal points are given as the two existing incinerators, the 
rd and the Pulaski. One additional alternative, that of rail- 


haul, ts also included. The unit cost used as based on the 
ly mentioned data accounts for a rail-haul of approximately 
Therefore, it is not necessary to locate the rail-haul disp 
facility. The x-y coordinates for the two existing facilit 
the unit cost of operation as given in the 1965 Public Work 
of the City of Baltimore are shown in Table 5-8. Refer to 
for location in relation to the collection area. 


Table 5-8. Characteristics of Baltimore Incinerato 


Name 


Location 

Capacity 
(tons /day) 

Reedbird 

- #3 

x * 26,000 

y = 45,000 

600 

Pulaski - 

#4 

x * 49,000 

y =* 26,000 

800 


From Anon. (1966a), Baltimore Public Works Department Annu 



D. Results of the Analysis 


1. Explanation of Symbols 

1 / 

A total of 60 runs were made with the facilities locatii 
model and the data case of section C. These runs represent a bi 
spectrum of changes in parameters and conditions in order to she 
the sensitivity of the system and to estimate the order of magn: 
>f the cost of proposed system policy and operational changes, 
runs are summarized in Table 5-9. The headings of each column s 
ixplained and referenced in the following manner: 

Type of Transfer Station - For most runs, seven transfer 
station site alternatives were allowed, A, B, C, D, E, F and RR. 
liese symbols and the locations of the sites have been shown 
►reviously in Table 5-6. It was assumed in most runs that the s 
apacity facility was proposed at each site. Any deviation from 
his policy, or in the number of sites available, is noted in th 
.EMARKS column beside the run. 

Collection Frequency - Two cases were considered: twice 
reek (2) and three times a week (3). 

Collection Costs - These are the costs of picking up the 
rolid waste within the collection areas. The basis for estimati 
:his cost is shown in section C (3) of this chapter. Units are 
lollars per week. 

-^Each run took an average of 45 seconds, or for the 60 runs a t< 
of 45 minutes of execution time on the IBM 7094, 


Transfer Cost from Collection - This represents t 
tation costs incurred by the collection vehicles in haul 
collected load from the area where it was picked up to t 
where it will be discharged, and the return trip to the 
area. The point of discharge in most cases may be eithe 
facility or a final disposal point. In most cases when 
stations were built not all collection vehicles used the 
still economically feasible for some collection subareas 
directly to final disposal. Units are dollars per week. 

Facility Costs - Cost per week of building and 0 { 
facility of a given capacity. Computations to arrive at 
figures are shown in Table 5-7. 

Transfer from Transfer Stations - This category i 
cost incurred using specialized tractor trailer transpoi 
to move wastes from the transfer site to the final dispc 
rail-haul is used, such transfer costs do not appear in 
category but as part of the disposal cost, since there j 
break apart the unit cost figure for rail-haul and dispc 
are dollars per week. 

Disposal Cos ts - This reflects the unit cost of c 
all cases where local disposal was chosen, the Pulaski i 
was chosen over the Reedbird incinerator because of its 
location with regard to the northwest quadrant that make 
study area, and because of the lower unit costs of dispc 
newer Pulaski incinerator. In the case where rail-haul 



Table 5-9. Description of Runs with Baltimoi 



Type of 

Collec- 

Collec- 

Transfer 

Fac'l’ty 

Trs 

Transfer 

tion 

tion 

from 

Costs 

i 

Station 

Freq. 

Cost 

Collec- 

Dollars 

Tre 

Proposed 

Times 

Dollars 

tion 

per Wk 

Sts 

Tons /Wk 

per Wk 

per Wk 

Dollars 


Do] 




per Wk 


. pe 




Load Load Total Total Cost Remarks 

Sent to Directly Load Dollars/Wk 

Fac'l't ’ s to Tons /Wk 

Tons/Wk Disposal 
Tons /Wk 


1060 

368 

1428 

15,599 



0 

1428 

1428 

17,399 



600 

828 

1428 

16,844 



900 

528 

1428 

16,606 



900 

528 

1428 

15,598 



600 

600 

228 

1428 

15,869 



1060 

368 

1428 

16,620 



1093 

335 

1428 

16,754 



1428 

0 

1428 

17,063 

Av'ge haul 
10 miles 

distance 

0 

1428 

1428 

18,419 

Av'ge haul 
10 miles 

distance 

0 

1428 

1428 

19,449 

Av'ge haul 
12 miles 

distance 

1428 

0 

1428 

17,349 

Av’ge haul 
12 miles 

distance 

1428 

0 

1428 

17,749 

Av'ge haul 
14 miles 

distance 

1428 

0 

1428 

19,855 

Av'ge haul 
14 miles 

distance 

1428 

0 

1428 

19,855 

Av'ge haul 
16 miles 

distance 

1428 

0 

1428 

18,149 

Av’ge haul 
16 miles 

distance 

0 

1428 

1428 

16,664 




Table 5-9. Description of Runs with Baltimore Data < 


Type of 

Collec- 

Collec- 

Transfer 

Fac'l'ty 

Trans 

Transfer 

tion 

tion 

from 

Costs 

frc 

Station 

Freq. 

Cost 

Collec- 

Dollars 

Trans 

Proposed 

Times 

Dollars 

tion 

per Wk 

Stat i 

Tons/Wk 

per Wk 

per Wk 

Dollars 


Dolls 




per Wk 


per 



Load 
>ent to 
r ac'l , t 1 a 
Fons/Wk 

Load 

Directly 

to 

Disposal 
Tons /Wlc 

Total 
Load 
Tons /Wk 

Total Cost 
Dollars/Wk 

Remarks 

1093 

335 

1428 

15,722 


1428 

0 

1428 

16,031 

Av'ge haul distance 
10 miles 

0 

1428 

1428 

17,815 

Av'ge haul distance 
10 miles 

1428 

0 

1428 

18,791 

Av’ge haul distance 
12 miles 

1428 

0 

1428 

16,481 

Av'ge haul distance 
12 miles 

1428 

0 

1428 

16,931 

Av’ge haul distance 
14 miles 

1428 

0 

1428 

18,791 

Av’ge haul distance 
14 miles 

1428 

0 

142 8 

18,791 

Av’ge haul distance 
16 miles 

1428 

0 

1428 

17,381 

Av’ge haul distance 
16 miles 

1202 

369 

1571 

17,592 

Increase waste load 
by 10% 

1315 

402 

1717 

19,134 

Increase waste load 
by 20% 

1427 

436 

1863 

20,670 

Increase waste load 
by 30% 

1534 

470 

2004 

22,149 

Increase waste load 
by 40% 

1649 

504 

2153 

23,717 

Increase waste load 
by 50% 


Table 5-9. Description of Runs with Baltimore Data i 


Type of 

Collec- 

Collec- 

Transfer 

Fac'l'ty 

Trans 

Transfer 

tion 

tion 

from 

Costs 

frc 

Station 

Freq. 

Cost 

Collec- 

Dollars 

Trans 

Proposed 

Times 

Dollars 

tion 

per Wk 

Statj 

Tons/Wk 

per Wk 

per Wk 

Dollars 


Do Ilf 




per Wk 


per 


'lty Load Load Total Total Cost Remarks 

It Sent to Directly Load Do liars /Wk 

Facfl*t 1 s to Tons /Wk 
Tons/Wk Disposal 
Tons /Wk 


1756 

538 

2294 

25,198 

Increase waste loac 
by 60% 

1800 

639 

2439 

26,737 

Increase waste loac 
by 70% 

1800 

706 

2586 

27,869 

Increase waste loac 
by 80% 

1800 

931 

2731 

29,818 

Increase waste load 
by 90% 

1112 

316 

1428 

15,659 

Transfer vehicle sp 
20 mph 

1135 

293 

1428 

15,323 

Transfer vehicle sp 
30 mph 

1147 

281 

1428 

15,209 

Transfer vehicle sp 
40 mph 

1093 

335 

1428 

16,604 

Vary collection rat 
by -10% 

1093 

335 

1428 

14,999 

Vary collection rat 
by +10% 

1093 

335 

1428 

14,349 

Vary collection rat' 
by +20% 

1500 

653 

2153 

22,270 

B is only alternate 
Increase waste lo« 
by 50% 

2100 

775 

2875 

24,516 

B is only alternati* 


Increase waste lo< 
by 100% 


Table 5-9 


Description of Runs with Baltimore Data 


Run 

No. 

Type of 
Trans fer 
Station 
Proposed 
Tons /Wk 

Collec- 
tion 
Freq. 
Times 
per Wk 

Collec- 
tion 
Cost 
Dollars 
per Wk 

Transfer 
from 
Collec- 
tion 
Dollars 
per Wk 

Fac'l'ty 
Costs 
Dollars 
per Wk 

Tran 

fr 

Tran 

Stat 

Doll 

per 

44 

1500 

3 

9,361 

1620 

1010 

7 

45 

1500 

2 

8,103 

1498 

2020 

8 

46 

1500 

3 

9,361 

1230 

2020 

5 

47 

1500 

2 

8,103 

1754 

1010 

9 

48 

1500 

3 

9,361 

1780 

1010 

7 

49 

1500 

2 

8,103 

1846 

1265 

7 


1500 

3 

9,361 

1620 

1265 

7 

51 

1500 

2 

8,103 

1846 

1518 

7 

52 

1500 

3 

9,361 

1620 

1518 

7 

53 

1500 

2 

8,103 

4564 

0 


54 

1500 

3 

9,361 

4040 

0 


55 

1500 

3 

9,828 

1700 

1010 

8 


Remarks 


y Load Load 

Sent to Directly 
FadlVs to 

Tons/Wk Disposal 
Tons /Wk 


Total 

Load 

Tons/Wk 


Total Cost 
Dollare/Wk 


1093 

335 

1428 

16,752 

B is only alternative 
Increase waste load 0% 

746 

411 

271 

1428 

16,471 

Political jurisdiction 
case 

700 

415 

313 

1428 

17,200 

Political jurisdiction 
case 

1044 

384 

1428 

15,805 

Deny B as an alterna- 
tive 

1428 

0 

1428 

16,863 

Deny B as an alterna- 
tive 

1093 

335 

1428 

15,977 

Increase fixed costa 
by 25% 

1093 

335 

1428 

17,009 

Increase fixed costs 
by 25% 

1093 

335 

1428 

16,230 

Increase fixed costs 
by 50% 

1093 

335 

1428 

17,262 

Increase fixed costa 
by 50% 

0 

1428 

1428 

16,664 

Increase fixed costs 
by 100% 

0 

1428 

1428 

17,399 

Increase fixed costs 
by 100% 

1148 

351 

1499 

17,539 

Increase waste load 
by 5% 

1202 

369 

1571 

18,321 

Increase waste load 


by 10% 


Table 5*9, Description of Runs with Baltimore Data t 


Type of 

Collec- 

Collec- 

Transfer 

Fac’l’ty 

Trans 

Transfer 

tion 

tion 

from 

Costs 

frc 

Station 

Freq. 

Cost 

Collec- 

Dollars 

Trans 

Proposed 

Times 

Dollars 

tion 

per Wk 

Statj 

Tons/Wk 

per Wk 

per Wit 

Dollars 


Doll* 




per Wk 


per 




ad Load Total Total Cost Remarks 

t to Directly Load Dollars/Wk 

L'fc's to Tons/Wk 
s/Wk Disposal 
Tons/Wk 


384 

1641 

19,103 

Increase 
by 15% 

waste 

load 

1499 

1499 

18,270 

Increase 
by 5% 

waste 

load 

1571 

1571 

19,132 

Increase 

waste 

load 


by 10% 

Increase waste load 
by 15% 


0 


1641 


1641 


19, 998 


as the disposal method, the disposal costs include the 
hauling to the distant disposal point. Where local dis 
chosen, the cost includes only the actual disposal proc 
sources of these costs are shown in section C (6) of th 
and are in dollars per week. 

Facilities Built - If transfer sites were chosen 
models, this indicates at which sites transfer stations 
given capacity would be built. 

Load Sent to Facilities - Tons of waste per week 
collection areas that would be sent to transfer facilit 
than directly to final disposal. 

Load Directly to Disposal - Tons per week sent d 
disposal without trans-shipment to a transfer site. 

Total Load - The total waste load in tons per we 
Represents the sum of the previous two categories. 

Total Cost - The sum of the collection, transfer 
collection, facility, transfer from facilities and disp 
Units are dollars per week. 

2. Explanation of Runs 

Estimation of Present System Costs and Verificat 
important aspect of the study is to establish some bene 
which the results of changing the system may be measure 
their efficiency. The present system runs with two tiir 
collection frequency so the computation of this cost wi 




bench mark but also as a check on the validity of the 
nee published figures are available for the cost of running 
1 system (Anon, 1966a). Run number 17 represents this 
ion since its data inputs are the best estimates of all 
lection frequency is twice a week and no transfer stations 
ad. The yearly cost of collection by the city from the 
division for the year 1965 was $722,000. The model 
lis cost as $12,666 per week, or $688,632 per year, which 
arence of 4.5 percent. The difference between the two 
is slight and may be explained by the fact that 1960 
\ estimates were used in the model while the actual cost 
presents the cost of collecting from the 1965 population 
le region. Since the northwest division of the city, 

•ly at its outermost fringes, is gaining in population, it 
ixpected that the model estimates would be slightly lower. 

imation of a Three Times a Week Collection Frequency Bench 
run number 2, an estimate was made of the cost of three 
ek collection frequency without transfer stations, with a 
total system cost of $17,399 per week. When compared to 
a week bench mark of $16,664 per week this shows only a 
it increase in cost to increase in collection frequency, 
his estimate is really only a lower bound on the actual 
tree times a week collection and thus should be viewed with 
•re are several assumptions buried in the calculation of 


the cost of three times a week collection that should be 


stated. The first is that the routes for the vehicles wi 
changed from those of the twice a week collection vehicle 
insure that the truck is close to capacity before it make 
to the transfer or disposal point. As noted, there is no 
mentation on a vehicle which tells the driver the present 
his vehicle. Routes are usually assigned in terms of def 
with the driver returning when he reaches a certain place 
in three times a week collection, less waste is generated 
area as there has been less time Bince the last collectio 
vehicle must travel farther to get a complete load. The 
vehicle is an important factor in the calculation. The a 
load for twice a week collection was 9000 pounds. Should 
average load for three times a week collection drop becau 
are improperly designed, say by 10 percent, the cost dif 
between twice and three times a week would rise to 8 perc 
the average collected weight drop by 25 percent, the cost 
between twice and three times a week collection would ris 
cent. Thus extreme care must be taken in designing the r 
make sure that enough load is picked up before transporta 
factor to be watched in the design of the routes is the t 
Most operations work under a M no overtime" constraint, wh 
that at a certain time the truck will return, regardless 
present load. Since it takes longer to collect under a t 
a week policy careful attention must be paid to the time 


a. Models as suggested in Chapter IV on routing would be 
ssign routes correctly. 

second assumption that has been made in these computations 
he waste load is the same for twice a week and three times 
llection. This may not be a valid assumption, but no 
evidence exists to disprove it. However, work by Quon et 
) based on actual observations in a controlled experiment 
5 shows significant increases in waste loads when col len- 
iency was increased from once to twice a week. They found 
i increases in the order of 30 to 50 percent. To test the 
waste load increase on the three times a week collection 
m, runs 58-, 59 and 60 were made. These runs represent 
‘s a week collection with no transfer, and a waste load 
if 5, 10 and 15 percent respectively. These runs show that 
ition of the difference between twice and three times a 
cfcion is quite sensitive to the assumption made about 
. These computations are shown in Table 5-10, Great care 
taken in estimating the difference between twice and three 
ek collection, as the cost difference is very sensitive 
i waste load estimates. 

rever, in measuring and comparing differences between indi- 
te when both have the same collection frequency, the type 
:ion made does not have a serious effect on the calculation 
same assumptions appear in both. The remainder of the runs 
a collection frequency of three times a week are all made 


with the assumption that the average truck load does not c 
from twice a week collection and that individual waste lo« 
generated does not change. 


Table 5-10. Different Waste Load Assumptions and Their Ei 
the Percent Difference Between Twice and Thrc 
a Week Collection with no Transfer 


Waste Load Assumption 

Total System Cost 
Dollars per Week 

Cost Diff* 
Between Twit 
and Three 1 
Week Coll« 

Same as 2 times a week 

17,399 

4.5/ 

5% greater 

18,270 

10.0? 

10% greater 

19,132 

15.0? 

15% greater 

19,998 

20.0? 


Is Transfer Feasible ? - With the best data availat 
various cost parameters and for the cost of facilities, a 
runs was carried out with different capacity transfer stal 
the proposed transfer sites to show if transfer stations ! 
cost of collection. In every case they did. The results 
analysis also give some interesting insights into the sysl 
3, 4, 7 and 8 were made for three times a week collection 
different sized transfer stations ranging from 40 percent 



; weekly waste loaf 1 (600 tons per week) to over 100 per- 
ons per week). Similarly, runs 6, 5, 1 and 18 were 
twice a week collection rate and the results of all 
re presented in Figure 5-4, First, notice that beyond 
per week alternative, the cost function for the higher 
tions is more favorable and quite flat, indicating little 
to size within that range. More important, for each 
y one transfer site was chosen and for all six cases 
s B. Thus this preliminary analysis would seem to favor 
esent system. In all cases, the building of only one 
tion indicated that some waste load was still going 
the final disposal point without transfer. Both the 
D ton capacity per week facilities were used at less 
/ with about 40 percent of the waste in both cases going 
disposal. Later runs will show that this proportion is 
/e to the location of the final disposal points. The 
>e times a week collection with transfer is comparable 
sent two times per week collection without transfer, 
it ions will show how comparable these costs are if 
iould be increased for the higher frequency case, 
le smaller capacity (600 tons per week) stations, a 
rernative (A) was chosen and in one case two alternatives 
;e chosen. In general, the difference between the solu- 
msfer stations and the one without is in the range of 


4 to 7 percent with the two times a week collection freq 
allowing a better decrease. The absolute difference bet 
transfer and non- transfer solution is a rough measure of 
facility costs can increase before transfer is no longer 
More will be said about this later. 

Effect of Increased Haul Distance to Final Dispos 
present average haul distance from the northwestern divi 
disposal at the Pulaski incinerator was shown to be 8 mi 
by Truitt (1968) and was verified by this author. Now a 
the haul distance was increased. Such a case is possibl 
the city expands, the tendency is to move waste disposal 
in order to avoid creating a nuisance to nearby neighbor 
because of the lack of availability of unused sites near 
Runs were made to see what the effect would be if the di 
points were moved so that average haul distance was incr 
12, 14 and 16 miles. In this analysis there is an addit 
possibility which was not considered in Truitt’s work, 
haul to far distant disposal points. As average haul di 
increases, transfer facilities become more and more favo 
finally the most dominant choice is rail-haul. A set of 
made to see how moving the disposaL point would affect a 
and without transfer stations. Runs 10, 11, 14 and 15 r 
case of increasing haul distance without allowing transf 
three times a week collection frequency. Runs 9, 12, 13 


tsfer under the same conditions. Figure 5-5 shows the 
: these runs compared to the cost of the rail-haul alter- 
,t is evident that as the haul distance increases, 
iecome8 more favorable to non-transfer. In fact with non- 
rail-haul becomes feasible at 12.9 miles average haul 
If transfer is allowed, rail-haul is not a suitable 
e until almost 26 miles average haul distance. This 
that as pressure mounts to move disposal out farther from 
a transfer facility system would preclude rail -haul for 
given that the relative costs of the alternatives would 
same. Further, the transfer site chosen in each of the 
he same site C, indicating that the site selection in 
i$ quite stable. Also, as soon as the disposal point has 
only 2 miles, all waste in the system goes through the 
acilities and none is taken directly to the disposal point, 
e reason C was chosen for this case while B was chosen for 
us cases with present location of facilities. The choice 
and C favors B only slightly, and when conditions are 
cause more load to go through a facility, C becomes the 
native. 

Figure 5-6, results of runs 19, 20, 21, 22, 23, 24, 25 and 
tted to show how a twice a week collection system would 
r the same haul distance changes. The results are quite 
Lth C being chosen as the transfer site and the rail-haul 
2 becoming feasible at a slightly smaller haul distance 


than in the previous case. 

Thus it would appear that the questions concerning transi 
and rail-haul are intimately linked. Under present conditions, 
rail-haul is not a feasible alternative. At a cost of $18,791 i 
week for the present system it represents a 12.8 percent increaj 
cost over the system with no transfer, and an 18.3 percent incr« 
over the system with transfer. However, if haul distance shoulc 
increase, rail-haul would become a feasible alternative in a fa: 
short time if no transfer is established. If transfer facility 
are built they reduce costs to the extent that rail-haul would i 
be feasible for a considerable time. As for transfer itself, it 
does exhibit some, but not extraordinary, savings over non-trans 
for the present system. For future conditions, however, these 
savings will increase, thus making transfer even more favorable, 

Sensitivity to other Parameters - Other runs were made tc 
show how sensitive the system was to estimates of the various 
parameters. Since waste load estimation in three times a week 
collection has been shown to have an effect on the solution, se\ 
additional runs were made to study this situation. Runs 55, 56 
57 represent runs with transfer facilities allowed, three times 
week collection frequency and an increase in the waste load of !. 
and 15 percent, while runs 58, 59 and 60 represent the system 
without transfer. The results are shown in Table 5-11. The cot 


of the transfer station decision is even more sensitive than th« 


e 5-11, Different Waste Load Assumptions and Their Effect on 
Percent Difference Between Twice and Three Times a 
Week Collection with Transfer 


ease in 
a Load 

11 

Increase in Cost 
with no Transfer 
(%) 

Increase in Cost 
with Transfer 
<%) 

Site 

Alternative 

Selected 

0 

4.5 

6.4 

B 

5 

10.0 

11.4 

B 

10 

15.0 

16.6 

B 

15 

20.0 

21.5 

B 


iransfer solution. However, in all cases the site selected was 
idieating that only the cost estimate but not the site chosen 
msitive to this parameter. Twice a week collection was tested 
.ts sensitivity to waste load increase in runs 27, 28, 29, 30, 
12, 33, 34 and 35, where the waste load was progressively 
;ased until it almost doubled. Results of these runs are shown 
.gure 5-7, Even though waste load was increasing, the solution 
responded linearly and the same alternative B, was chosen for 
solutions . Even when waste load was almost double the transfer 
Lon capacity it was not optimal to build a second facility at 
Lferent site. 


Runs 36, 37 and 38 show how sensitive the solution is to the 



speed of the transfer vehicle from the transfer facility to t 
disposal point. Results are plotted in Figure 5-8. Present 
mates place this speed at 16 miles per hour since this is the 
approximate speed in traffic of the collection vehicles, and 
transfer vehicles are larger tractor trailers which must opet 
under largely the same traffic conditions. The plot indicate 
an increase in the traffic speed from 16 to 30 miles per hout 
bring about a decrease of about 7.5 percent in the cost of tV 
solution. However, it would be expected that the cost of the 
remedial measures necessary to bring about this increase in £ 
would more than outweigh any savings in transportation costs. 

An attempt was also made to see how collection rates £ 
the solution. In runs 39, 40 and 41 collection rates are vat 
by -10 percent, +10 percent and +20 percent. Increasing or c 
the collection rate by 10 percent leads to about a 5 percent 
in the solution cost, which falls into the range of being son 
sensitive. Thus care should be taken in estimating these rat 
Further, work rules that tend to increase collection rate wit 
too great an increase in cost or inexpensive time saving devd 
would also bear some investigation. 

The effect of changing the estimates of the facilities 
is also of interest. This gives some feeling for how careful 
estimates have to be made for preliminary studies. Runs 49, 
52, 53 and 54 show how fixed cost increases of 25, 50 and 10C 
would affect the solution. At about a 75 percent increase ii 


costa, transfer facilities are no longer feasible, but up to that 
point the same facility was chosen each time. A rough estimate of 
how much facilities costs could increase could have been found frc 
earlier runs simply by looking at the absolute difference in cost 
between the transfer and non-transfer solutions. If the facility 
cost increase exceeded this difference, transfer was no longer 
feasible. 

Political, Aesthetic and Regional Constraints - Many times 
constraints occur which preclude the use of an alternative that 
would otherwise seem to be the best from an economic point of view 
Examples are political constraints when a site might better be put 
to use for other municipal functions, aesthetic constraints where 
the neighborhood around a site strongly opposes a facility there, 
and regional constraints which do not allow for cooperation betwee 
different political subdivisions. Consider first an example where 
a particular site is not allowed to be used. The results of the 
analysis to this point show that B is the best transfer site. Run 
47 and 48 show what happens when this site is excluded from the se 
of alternatives. In the case of two times a week collection, Site 
is chosen and the increase in the cost 1 b about $60 per week. For 
three times a week collection the site chosen is site C and the 
difference is about $110 per week. Both of these changes are 
minuscule compared to the total cost of operating the system. Thi 
indicates that there may be readily available alternatives to the 


best site that will not raise the cost of the solution signifi 
Thus the cost of satisfying additional constraints which prohi 
certain sites is not particularly high, and the decision make* 
a great deal of leeway in choosing between alternatives when c 
criteria must also be considered. 

The effect of regionalization will only be treated brief] 
this example, as neither the time nor data were available for 
assembling an extremely large regional case. However, the quc 
can be looked at implicitly by assuming that the northwestern 
division is a unified region, and asking how much costs would 
increase if it were subdivided into two non-cooperating regiot 
Suppose the region were divided into two so that the northern 
could use only transfer sites A, D and E, and the southern are 
could use only sites B, C and F. The results of these runs at 
shown in Table 5-12. In this case subdividing the region sti3 
makes transfer feasible, but only at an extremely slight advai 
over no transfer at all. This would suggest that regionalizat 
would mean increased savings. Although this analysis was not 
extended to the entire city, it might be expected that each oi 
administrative divisions would not have a transfer facility, t 
only two or three would be built. 
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Table 5-12. Effects of Regional Cooperation 


u u 

•H <i> 
p iw 
•H CO 

n e 

n} 

c u 
u e-< 

V 


CO 3 

<y o 
? ^ 
►O 4J 
+J -H 

n 

O 

55 


M 

*J QJ 

o m 


4-1 

3 

O 

.c 


C ^ 

o u 

•H -H 
CO J? 
•H 

> « 

•H 4-1 

O O 


U 
a i 

C l-i 4-i 

M +J to 

<u w C 

4-« *H fl) 

S Q 6 
S o 


*<5 


4J 

u 

O CO 

55 *H 


g * 

■H O 
4-1 C 


3 

f~4 D* 
r-4 <U 

o u 

o tM 





<u 

<r 

Ov 

o 

vO 

O' 


vO 

co 

\ 

a 

— 

u 

vO 

r- 

CO 

T-H 

f-H 

o 

co- 


a 



<n 



■w 




ca 

CQ 

W 



^5 



0) 

o 

CM 

<D 

CM 

m 

£ 

r-. 

** 

c-» 

A 

*J 

in 

vO 

CO 

*-< 

i — 1 

o 

</> 


a 




Q 

f*4 

CO 



01 

•o 

T3 

4-1 

C 

C 

•r4 

al 

nj 

CO 




« 

o 




0 

1—4 

o 

a) 


o 

£ 

a 

CM 
* 1 

•M 

vO 

I s * 

<0 

H 

H 

o 

</> 


O 




CM 

CO 



Summary of the Experimental Analysis 


The following general statements about the system can be 
i, based on the rough analysis with estimated data for the 
:hwestern division. These runs have been intended as a demon- 
ition of the techniques of this thesis, rather than a detailed 
Lysis of the system. The actual findings would need verifica- 
1 with better cost data before they could be viewed as 
>letely reliable. 

1. The building of transfer facilities appears to be feasible 
and would result in an annual savings, under present condi- 
tions, of 7 percent in the total cost of operating the solid 
waste collection system. This saving would be expected to 
increase in the future, as the system expands and grows, 

2. For the northwestern division only one transfer facility 
would be built, with site B chosen as the location. Site C 
is also a good alternative and selecting C over B would 
show little disadvantage. 

3. The most favorable size of the facility to be built is 
1500 tons per week capacity. At this size some of the 
waste material will be taken to transfer and some will 
still be taken directly to the disposal point. These 
assignments are shown by the model. 



A rail-haul alternative would cost about 12 percent more to 
operate than the present system and about 18 percent more 
than the system with a transfer facility. Haul distance 
would have to increase an average of from 3.5 to 5 miles 
for rail -haul to be favorable without transfer, and from 
14 to 18 miles to be favorable with transfer. 

The cos t of changing from two times a week frequency of 
collection to three times a week is difficult to estimate 
because of the sensitivity of this calculation to estimates 
of how much the waste load would be with three times a week 
collection. Estimates of the added cost range from 4.5 per- 
cent if there is no increase in the waste load, up to 
21.5 percent with a 15 percent increase in the waste load. 
All of these figures are based on the assumptions that 
routes for vehicles under three times a week collection be 
redesigned for efficiency purposes. It is strongly advised 
that more studies of the waste load question in this region 
be carried out before any firm decision be made on changing 
collection frequency. 

While differences in assumptions make estimating the cost of 
some programs difficult, the site selection B, for a transfe 
facility is remarkably insensitive to parameter change. Re- 
gardless of the collection frequency, site B is the best sit 
for the transfer facility. 


7. The system shows some sensitivity to collection rates, 
which indicates that care should be taken in measuring 
them and investigation of means of altering them will be 
of some help. 

8. If for some reason site B could not be used, site G first 
and then site A are good alternatives, and it would cost 
little extra to change the selection to them. The effects 
of regionalization are such that the northwestern division 
is better operated as one division than as two. It is 
further expected that extension of this type of analysis 
to the entire city might show even greater efficiency with 
transfer, and this is recommended as the next step in the 
analysis procedure. 

9. The amount of computer time necessary for this analysis of 
60 runs was 45 minutes, using the IBM 7094. At commercial 
rates of $500 per hour, this represents a cost of 
approximately $375. 



CHAPTER VI. CONCLUSIONS AND EXTENSIONS 


The goal of this thesis has been the development of too 
gain a better understanding of a large-scale public system. T 
it is on the basis of the results of the example rough analysi 
the Baltimore, Maryland, solid waste collection system in Chap 
that the success of such an undertaking should be measured. W 
most of the data used for the analysis were less than elegant 1 
obtained, a firm body of facts and conclusions emerge from the 
exercise. Answers to many of the basic questions with which t 
analyst is concerned were found. Most important, a feeling fo 
sensitivity of the system to many of its parameters was establ 
The decision to change from twice a week to three times a week 
collection was shown to be quite sensitive to estimates of was 
generation under the two regimes. Thus, it is recommended tha 
additional study in this area be carried out before a decision 
made about a policy change. On the whole, however, the proble 
choosing transfer facility location appears to be quite stable 
the same location being elected under a severe change in syste 
conditions . 

Some drawbacks in the analysis were encountered because 
the need to subdivide the problem in order to solve it. At or 
point it was necessary in the study of facility location to as 
that the new routes for vehicles would be set up efficiently i 
service policy whould change. However, the model could not 


determine how this was to be done, and the problem would have to 
looked at as a separate vehicle scheduling problem as described 
Chapter XV. 

From a computational point of view the analysis in Chapte 
was a success. Sixty different problems, each involving general 
data for forty waste sources and considering nine intermediate 
alternatives, were solved in about three-quarters of an hour on 
IBM 7094. This represents a cost at commercial rates of less th 
400 dollars, or about two weeks' salary for a medium level 
engineering analyst. In terms of the information generated, it 
seems like a considerable bargain. 

Looking at the individual models developed, the large-sca 
’’flow of goods” models developed for the facility location probl 
of Chapter II and the multi -commodity truck assignment problem o 
Chapter III proved to be the most successful. When programmed f 
the computer, they were of sufficient size and speed to handle a 
very large complex system and find solutions quickly. Since muc 
of the value of the models is in repeated solution with changing 
parameters to determine sensitivity, efficient solution time is 
great importance. The vehicle scheduling models of Chapter IV d 
not display the same qualities. People working in this field ha 
uniformly found that securing an optimal solution to even a smal 
scale problem in a reasonable amount of time was difficult to 
impossible. Although the author's work in a problem area of thi 
general category extended theoretical concepts involving the use 


n one terminal, the same type of problem occurred. Thus, 
e no existing optimization techniques that will allow the 
ation of a city-sized routing problem. The development 
Stic methods has advanced, however, and feasible solutions 
e complex problems can be found quickly. No one has yet 
ated the difference between a good feasible solution and 
al solution in vehicle routing to determine if the search 
mal procedures is worthwhile. Meanwhile, sensitivity 
using heuristic rather than optimal procedures must be 
ith extreme caution and care. 

arhaps the problems of routing would benefit from a shift 
sis from discrete problems to continuous problems. Vehicls 
ng is viewed as a problem of finding routes between 
points or in terms of a network example, in finding which 
travel to visit all nodes. The rapidly increasing 
jrial nature of the problem is what makes computation of 
solutions so difficult. Studying public services in city 
itworks would reveal a continuous distribution of demand 
arc with a requirement that all arcs be traveled. For 
iases, this turns out to be a problem that is remarkably 
quick to solve. The Chinese postman problem asks how a 
is route for one vehicle can be found through a network 
'els all arcs and minimizes total distance traveled. The 
l of these techniques to problems of more than one service 


vehicle, each with a capacity constraint, while retaining the 
property of efficient solution would mean a dramatic increase in 
the type of analysis done on vehicle routing in solid waste 
collection. 

The overall conclusion of this thesis is that models such 
as developed, limited though they may be by applying only to sub- 
problems of the total collection system, can still provide a 
great deal of insight into the system. The types of questions 
amenable to study using these techniques, provide a wealth of 
information to the decision maker charged with the operation of 
the system. It is hoped that such models will become increasingl 
valuable tools in the management and planning of large-scale soli 
waste collection systems. 
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